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ABSTRACT 


Work  is  continuing  primarily  in  gas  phase  turbulent  mixing  and 
chemical  reactions.  The  liquid  phase  work  to  date  is  in  its  final 
stages  of  being  analyzed  and  documented  for  dissemination  in  the  form  of 
archival  publications.  In  the  gas  phase  shear  layer  work,  our 
investigations  are  concentrating  on  shear  layer  free  stream  density 
ratio  effects,  finite  kinetic  rate  (Damkohler  number)  effects,  and  a 
design  effort  in  support  of  the  planned  extension  of  the  work  to 
supersonic  flows.  In  jet  flows,  progress  has  been  made  in  the  gas  phase 
laser  Rayleigh  scattering  techniques  developed  for  conserved  scalar 
measurements  down  to  diffusion  space  and  time  scales.  A  new  technique 
has  been  developed  under  joint  support  with  the  Ca3  Research  Institute 
that  permits  the  imaging  of  soc;  sheets  in  turbulent  flames  and  is  being 
used  to  describe  the  combustion  flame  sheets  in  methane  flames. 
Theoretical  work  in  progress  is  addressing  the  finite  chemical  rate 
problem  as  well  as  the  diffusion-limited  shear  layer  mixing  problem. 
Advances  in  our  data  acquisition  capabilities  during  the  last  year  are 
permitting  higher  temporal  resolution  measurements  to  be  taken  with 
digital  image  arrays. 
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1 .0  INTRODUCTION 

Work  is  continuing  primarily  in  gas  phase  turbulent  mixing, 
chemical  reactions  and  combustion,  while  the  work  in  liquid  phase 
systems  is  in  its  final  stages  of  being  analyzed  and  published.  A 
review  paper  has  just  appeared  (Broadwell  &  Dimotakis  1986)  summarizing 
these  early  results  and  their  implications  for  modeling  chemically 
reacting  flows.  Both  experimental  and  analytical/modeling 
investigations  are  in  progress,  as  well  as  an  accompanying  effort  to 
develop  and  augment  diagnostic  techniques,  on  an  as  needed  basis,  in 
support  of  the  experimental  program.  Progress  in  these  areas,  realized 
during  the  reporting  period  ending  15-Apr-86  will  be  described  below. 

2.0  HF  COMBUSTING  SHEAR  LAYER 

2.1  Heat  release  effects 

The  experimental  investigation  of  the  effects  of  heat  release  (and 
moderate  pressure  gradients)  on  equal  free  stream  density  shear  layers 
has  been  completed.  See  discussion  in  our  Interim  Report  for  the  period 
ending  15~Apr-1985  (Liepmann,  Broadwell  &  Dimotakis  1985,  section  2.4 
and  appendix).  A  documentation  of  that  work  has  appeared  during  the 
period  covered  by  this  report  in  the  form  of  a  GALCIT  Ph.D.  thesis 
(Hermanson  1985).  Minor  revisions  of  the  first  announcement  of  that 
research  (Hermanson,  Mungal  &  Dimotakis  1985)  are  in  progress  for 
publication  in  the  AIAA  Journal. 

2 . 2  Stability  analysis 

A  study  was  carried  out  to  determine  the  effect  of  large  density 
ratios  on  the  initial  stability  characteristics  of  the  plane  mixing 
layer.  In  this  study,  the  inviscid,  linear,  parallel-flow  stability 
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analysis  of  spatially  growing  disturbances  was  utilized  to  calculate  the 
range  of  unstable  frequencies  and  wave-numbers  in  a  non-uniform  density 
plane  shear  layer  taking  into  account  the  wake  component  of  the  initial 
velocity  profile.  This  work  has  been  submitted  for  the  AIAA  25th 
Aerospace  Sciences  Meeting  in  Reno,  Nevada,  January  1987  (see  appendix 
A).  The  most  important  result  is  that  when  the  high  density  >s  on  the 
low-speed  side,  there  exist  two  instability  branches  of  comparable 
growth  rate.  One  of  these,  the  shear  layer  mode,  leads  to  the  usual 
Kelvin-Helmholtz  roll-up  patterns.  The  second  branch,  the  wake  mode, 
corresponds  to  roll-up  patterns  that  resemble  those  found  in  pure  wake 
flows.  These  results  indicate  that,  under  the  right  conditions,  a  shear 
layer  of  non-uniform  density  may  transition  to  turbulence  more  like  a 
wake  than  like  a  shear  layer  through  the  Kelvin-Helmholtz  vortices. 
That  this  can,  in  fact,  happen  was  confirmed  in  an  acoustically-forced 
shear  layer  between  a  high-speed  stream  of  helium  and  low-speed  stream 
of  argon.  Both  shear  layer  and  wake  type  flow  patterns  were  generated 
depending  on  the  forcing  frequency.  See  figure  4  in  appendix  A. 


2.3  Damkohler  number  effects 

A  property  of  the  hydrogen-fluorine  reaction  system  that  has  proved 
to  be  most  useful  is  that  its  overall  rate  is  determined  by  the  hydrogen 
atom  concentration,  which,  in  turn,  can  be  controlled  by  the  initial 
nitric  oxide  concentration.  This  system  is,  therefore,  almost  ideal  for 
studying  finite  rate  chemistry:  the  potential  heat  release  can  be  small 
and  fixed  as  the  rate  is  changed  from  zero  to  values  so  high  that  the 
reaction  is  mixing  limited.  Previous  extensive  studies  (Mungal  & 
Dimotakis  1984,  Mungal,  Hermanson  &  Dimotakis  1985)  at  the  high  rate 
provide  a  set  of  results  against  which  those  obtained  at  "finite  rates" 
can  be  compared. 


In  the  past  year,  a  series  of  experiments  were  carried  out  in  the 
H2“f2  shear  layer  combustion  facility  using  this  reaction  to  investigate 
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the  effects  of  reaction  rate  in  a  high  Reynolds  number  shear  layer.  In 
addition,  the  finite  reaction  rate  was  incorporated  into  the  model 
proposed  by  Broadwell  and  Breidenthal  (1982)  and  which  has  been  used  in 
the  interpretation  of  the  earlier  H2-F2  high  reaction  rate  experiments. 
These  studies  are  described  in  Mungal  &  Frieler  (1985)  and  Broadwell  & 
Mungal  (1986),  which  are  attached  as  Appendices  B  and  C. 

The  primary  experimental  results  were  sets  of  temperature-time 
traces  recorded  at  eight  stations  distributed  laterally  across  the  layer 
at  one  axial  location.  From  these,  mean  temperature  profiles  can  be 
constructed  and  the  total  amount  of  product  in  the  layer  calculated.  In 
the  model,  the  parameter  which  arises  naturally  to  characterize  the 
reaction  rate  is  the  first  Dahkohler  number,  Dax>  the  ratio  of  two 
times:  the  mean  flow  time  between  the  measuring  station  and  the 
transition  location  and  a  characterisitc  reaction  time.  The  first  time 
is  2  x/ ( U 1  +  U2),  where  the  U's  are  the  free  stream  velocities  and  x  the 
distance  just  mentioned.  The  reaction  time  was  determined  from  the 
CHEMKIN  code  for  the  various  free  stream  reactant  concentrations. 

The  experimental  and  analytical  results  are  shown  in  figure  1  taken 
from  Broadwell  &  Mungal  (1986).  We  see  that,  for  the  fastest  rates,  the 
experiments  and  the  model  agree  in  showing  that  the  reaction  is  nearly 
independent  of  Dax,  i.e.,  nearly  mixing-limited  at  the  highest  rates. 
Finite  rate  effects  become  significant  for  Dax  -  30  (a  result  of 
critical  importance  for  the  design  of  the  supersonic  facility). 

In  addition,  the  changes  in  the  shapes  of  the  time  traces  and  of 
the  mean  profiles  with  reaction  rate  in  the  experiment  are  important 
clues  to  the  nature  of  the  turbulent  mixing  process.  The  results  of  a 
preliminary  evaluation  in  Broadwell  &  Mungal  (1986)  shows,  again,  the 
central  role  of  the  large  scale  structures  in  the  first  stages  of  the 
mixing  process.  The  observations  are  also  consistent  with  the  ideas 
underlying  the  model.  Interestingly,  these  ideas  explain  why  the 
time-traces  and  profile  shapes  for  a  slow  reaction  in  a  gas  are  similar 
to  those  for  a  fast  reaction  in  water. 
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With  regard  to  the  model  we  note  that  there  are  two  empirical 
constants  (one  of  which,  however,  can  be  estimated  from  the  turbulence 
theory  on  which  the  model  is  based).  For  the  same  values  of  the 
constants,  the  model  satisfactorily  accounts  for  the  effects  of: 

1 .  stoichiometric  mixture  (equivalence)  ratio  (Mungal  &  Dimotakis 
1984), 

2.  Schmidt  number  (Koochesfahani  1984), 

3.  Reynolds  number  (Mungal,  Hermanson  &  Dimotakis  1985)  —  not  so 
satisfactorily, 

4.  the  "flip"  in  liquids  (Koochesfahani  1984,  Koochesfahani, 
Dimotakis  &  Broadwell  1985), 

and ,  now , 

5.  the  DamkcJhler  number. 

When  the  same  ideas  are  used  in  a  model  of  fuel  jets  burning  in 
air,  it  explains  (Broadwell  1982): 

1  .  the  independence  of  flame  length  from  Reynolds  and  Schmidt 
numbers  at  high  Reynolds  number, 

2.  the  linear  dependence  of  flame  length  on  the  stoichiometric 
mixture  ratio, 


3.  the  dependence  of  nitric  oxide  generation  in  hydrogen  and 
methane  flames  on  Reynolds  number. 

Finally,  if  the  influence  of  strain  rate  on  the  chemical  reaction  is 
introduced , 
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4.  the  parameters  governing  blow-out  of  hydrocarbon  flames  can  be 
deduced  (Broadwell,  Dahm  and  Mungal  1984). 


These  results  suggest  that  the  single  idea  on  which  the  model  is 
based,  namely  that  molecular  mixing  takes  place  in  a  succession  of  a 
Kolmogorov-like  cascades  between  the  largest  and  smallest  scales  in  the 
flow,  embodies  most  (or  all)  of  the  essential  features  of  the  turbulent 
mixing  process.  The  model,  as  presently  formulated,  describes  the  mean 
behavior,  the  average  over  the  cascade  processes.  The  underlying 
processes  are  unsteady  and  the  observed  large  scale  flame  length 
fluctuations  are  the  best  evidence  of  the  individual  cascade.  A  more 
detailed  examination  of  the  results  of  several  experimental 
investigations  will  be  undertaken  to  determine  to  what  degree  this 
conclusion  is  justified. 


2.4  Chemical  kinetics  calculations  with  mass  addition 

In  our  extension  of  the  work  to  date  to  supersonic  flows,  we  would 
like  to  retain  the  considerable  advantage  of  being  in  a  position  to 
measure  the  extent  of  turbulent  mixing  to  a  molecular  scale  using  fast, 
irreversible  chemical  reactions  Accordingly,  the  need  for  fast  chemical 
kinetics  in  the  supersonic  shear  layer  flow  was  recognized  as  being  the 
primary  design  driver  since,  for  supersonic  flow,  the  residence  time  of 
the  reactant  gases  in  the  test  section  can  be  very  short  (order  of  1  to 
10  milliseconds).  It  is  this  consideration  which  motivated  our 
selection  of  the  F2/H2/NO  reaction  system  and  our  intention  to  have  the 
flow  test  section  as  near  as  feasible  to  atmospheric  pressure. 

Our  goal  is  to  be  able  to  perform  experiments  in  the  mixing-limited 
regime,  in  which  the  chemical  kinetic  rates  are  large  relative  to  the 
mixing  rate,  that  is  in  the  limit  of  infinitely  fast  chemical  kinetics. 
In  order  to  predict  whether  or  not  this  will  be  possible,  we  have 
resorted  to  a  two-step  analysis.  First,  we  have  assumed  that  the 
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experimental  results  of  Mungal  and  Frieler  (1985)  can  be  reasonably 
extended  to  supersonic  shear  layers.  This  study  indicated  that  the 
mixing  limited  regime  is  attained  in  flows  for  which  Dax  (Damkohler 
number  based  on  length)  is  greater  than  40.  Second,  we  performed 
computer  calculations  to  estimate  characteristic  chemical  times  for 
reactant  compositions  that  we  expect  in  our  supersonic  experiments. 
These  calculations  enable  us  to  estimate  the  maximum  attainable 
Damkohler  number  under  the  least  favorable  flow  conditions,  namely  those 
at  the  highest  flow  Mach  number  of  3.0  . 

Our  original  computer  calculations  followed  those  of  Mungal  and 
Frieler  in  that  they  simulated  a  constant  mass,  homogeneous  chemical 
reactor.  As  was  previously  reported,  those  results  suggested  that  our 
new  supersonic  facility  would  be  capable  of  a  mixing-limited  chemical 
kinetic  environment,  with  approximately  an  order  of  magnitude  safety 
margin.  We  are  presently  exploring  a  more  realistic  chemical  kinetics 
reactor  model  in  which  we  are  simulating  the  continuous  mass  (reactant) 
addition  that  results  from  the  entrainment  and  mixing  processes  and 
which  proceeds  concurrently  with  the  chemical  reaction.  The  new  results 
demonstrate  that  characteristic  chemical  times  can  be  significantly 
longer  when  mass  addition  is  included.  See  figure  2  .  These 
preliminary  results,  however,  indicate  that  it  may  still  be  possible  to 
operate  in  the  mixing-limited  regime,  albeit  with  a  somewhat  reduced 
margin  with  respect  to  the  fast  chemistry  requirement. 

There  are  two  caveats  that  must  be  stated  with  regards  to  the  fast 
chemistry  criterion  as  it  relates  to  the  supersonic  shear  layer.  First, 
it  is  not  clear  that  the  subsonic  results  of  Mungal  and  Frieler  can  be 
applied  directly  to  supersonic  shear  layer  flows.  The  turbulent 
interface  generation  under  supersonic  flow  conditions  could  be 
qualitatively  and  quantitatively  different  and,  as  a  consequence,  our 
estimated  margins  for  fast  chemistry  could  prove  optimistic.  Second, 
fast  chemical  kinetics  can  only  be  attained  with  relatively  high 
concentrations  of  reactant  gases;  for  example,  the  fluorine  gas 
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concentration  must  be  10  to  15  molar  percent.  This  results  in  elevated 
gas  temperatures  (in  excess  of  1300  K)  which  may  make  it  difficult  to 
study  the  mixing-limited,  supersonic  shear  layer  as  decoupled  from  the 
effects  of  heat  release.  These  are  important  issues  which  will  be 
addressed  in  the  course  of  the  research  that  will  be  undertaken.  They 
have  prompted  us  to  consider,  however,  the  additional  flexibility  that 
adding  supersonic  stream  gas  pre-heating  capability  to  the  design 
(subject  to  the  availability  of  funding  support)  will  permit. 


2.5  Supersonic  shear  layer  design 

The  GALCIT  Reacting  Shear  Layer  Combustion  facility  is  in  the 
process  of  being  re-designed  for  supersonic  flows.  The  new  facility  is 
being  designed  to  produce  a  bisonic  (supersonic/subsonic)  reacting  shear 
layer  capable  of  handling  the  reactants  of  the  F2/h2/N0  system.  The 
fundamental  design  considerations  behind  this  facility  were  described  in 
some  detail  in  our  last  report.  For  convenience,  the  supersonic  stream 
plumbing  schematic  (figure  3)  and  the  table  of  specifications  (Table  1) 
are  reproduced  here.  It  should  be  recalled  that  the  low  speed  (F2) 
plumbing  hardware  in  the  existing  laboratory  will  be  used  without 
significant  modification  in  the  new  facility.  Therefore,  our  design 
efforts  have  primarily  dealt  with  the  new  supersonic  leg  and  test 
section . 

Four  segments  of  the  facility  for  the  new  supersoinc  leg  have 
received  attention  since  our  last  report.  They  are:  the  reactant  tank, 
the  (possible)  ore-heater,  the  fast  acting  valve,  and  the  main  pressure 
regulator.  Although  none  of  these  items  are  ready  for  fabrication  or 
purchase,  it  is  hoped  that  only  one  more  design  iteration  for  each  will 
be  required. 

The  desirability  of  a  gas  pre-heater  was  explained  in  the  previous 
section.  The  current  design  calls  for  the  reactant  tank  to  be  wrapped 
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with  electrical  heating  tape  and  thermal  insulation.  The  total  heating 
capability  will  be  10  kW.  Prior  to  the  experiment,  the  tank  and  its 
packing  will  be  heated  to  maximum  temperature  of  600  K.  After  the 
desired  temperature  is  reached,  the  supersonic  gas  nr  'ture  will  be  added 
to  the  tank.  During  the  tank  discharge,  the  tank  ?nd  its  packing  should 
limit  the  temperature  drop  to  approximately  20  K. 

A  new  vendor  for  the  fast  acting  valve  has  been  located  (Fox  Valve 
Development  Corporation,  East  Hanover,  NJ ) .  This  valve,  which  wa3 
recommended  to  us  by  Paul  Waltrup  of  the  Applied  Physics  Laboratory  of 
Johns  Hopkins  University,  is  of  the  venturi-type  with  much  improved 
performance  over  the  original  device.  Opening  times  on  the  order  of  100 
ms  or  less  are  claimed  for  this  valve. 

Another  product  from  the  Fox  Corporation  could  potentially  serve  as 
a  main  pressure  regulator.  The  device  is  an  adjustable  area  venturi 
valve  which  could  be  combined  with  an  electronic  controller  so  as  to 
provide  mass  flow  regulation  and  hence  pressure  regulation.  The 
feasibility  of  such  an  actively  controlled  device  is  in  the  process  of 
being  evaluated.  After  the  feasibility  study,  the  relative  merits  of 
that  device  will  be  compared  with  those  of  the  commercial  pressure 
regulators  already  located,  and  a  final  device  selected. 


2.6  Density  Ratio  Effects 


By  selecting  the  free-stream  diluents  in  the  HF  shear  layer 
combustion  facility  it  is  possible  to  study  the  effects  of  density 
differences  on  entrainment,  mixing  and  reaction.  Keeping  the  velocity 
ratio,  flame  temperature  and  molar  specific  heats  constant,  experiments 
will  eventually  be  performed  for  the  range  of  density  ratios  1/8  <  s  <  8 
(s  =  P  2 / P 1 ) -  Pressure  gradients  due  to  displacement  effects  can  be 
removed  by  properly  adjusting  the  test  section  side  walls. 
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Several  important  questions  are  addressed  in  this  work.  Among 
these  is  the  dependence  of  the  shear  layer  volumetric  entrainment  ratio 
and  mixed  fluid  composition  on  density  ratio.  The  flux  of  fluid  into 
the  turbulent  region  can  be  estimated  from  the  orientation  of  the 
free-stream  streamlines  and  the  boundaries  of  the  mixing  layer 
(Dimotakis  1984,  Hermanson  1985).  The  ratio  of  the  fluid  flux  from  the 
two  streams  is  the  volume  entrainment  ratio  (Ev)  and  has  been  shown  to 
be  sensitively  dependent  upon  the  density  ratio  (Konrad  1976,  Mungal  & 
Dimotakis  1984).  Shown  in  figure  4  are  preliminary  data  for  Ey  for  a 
number  of  density  ratios,  along  with  a  similarly  derived  empirical  value 
from  the  Brown  &  Roshko  (1974)  data.  It  should  be  noted  that 
theoretical  analyses  by  Brown  (1974)  and  Dimotakis  (1984)  predict  a  /s 
dependence  on  density  ratio  for  this  quantity.  The  dashed  line  in 
figure  4  is  a  plot  of  the  theoretical  expression  (Dimotakis  1984) 

EV  =  +  0.68  ) 

where  r  -  U2/U1  is  the  velocity  ratio  which  was  set  in  the  range  of 
0.38-0.40  for  the  data  in  figure  4.  (note  that  Ey  >  1  for  p2/pi  =  1  ) . 
While  this  estimate  is  in  good  agreement  with  the  data  for  density 
ratios  s  >  1,  the  behavior  for  density  ratios  s  <  1  appears  to  be  less 
well  described  by  this  theory. 

The  anomalous  point  for  s  =  10  is  worth  elaborating  upon.  Detailed 
elsewhere  within  this  report  is  recent  work  by  Koochesfahani  and  Frieler 
which  demonstrates  the  existence  of  a  second  fundamental  mode  of 
instability  in  the  experimentally  realizable  shear  layer.  This 
corresponds  to  an  instability  of  the  splitter  plate  wake  velocity 
profile  arising  from  the  merging  of  the  boundary  layers  shed  from  the 
trailing  edge  of  the  splitter  tip.  The  predominance  of  this  mode  in  the 
cases  with  s  »  1  appears  to  be  responsible  for  significant  changes  in 
the  structure  and  dynamics  of  the  initial  roll-up  of  the  mixing  layer, 
which  seem  to  extend  for  a  significant  distance  downstream.  We  believe 
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that  factor  of  two  to  three  decrease  in  Ey  at  s  =  10,  over  the  value 
that  would  have  been  anticipated  by  any  reasonable  extrapolation  of  the 
rest  of  the  data,  is  a  manifestation  of  this  qualitative  difference  in 
shear  layer  behavior.  This  is  also  reflected  in  the  layer  growth  rate 
and  mean  stagnation  pressure  profiles,  which  also  indicate  that  the 
effects  of  this  mode  last  for  the  extent  of  the  HF  facility  test 
section.  Further  work  on  these  extreme  density  ranges  will  rely  on  our 
ability  to  manipulate  initial  conditions  to  control  this  instability. 

Using  chemical  reactions  at  stoichiometric  mixture  ratios  both 
large  and  small  compared  to  unity,  e.g.  <|>  *  8  and  $  =■  1/8  ,  it  is 
possible  to  estimate  the  volume  ratio  comprising  the  mean  mixed  fluid 
composition  (see  Koochesfahani  1984),  denoted  here  by  EvM  .  in  this 
work  this  is  derived  from  pairs  of  "flip"  experiments  performed  at  the 
same  density  ratio  at  high  and  low  free-stream  reactant  stoichiometric 
mixture  (concentration)  ratios  (<)>  -  8,  1/8).  Also  plotted  in  figure  4 
are  the  estimates  of  the  entrainment  ratio  derived  in  this  manner  (EyM) 
from  these  "flip"  experiments,  along  with  the  equal  density  experiments 
of  Mungal  &  Dimotakis  (1984,  gas  phase)  and  Koochesfahani  (1984, 
liquid).  Results  from  two  (non-reacting)  gas  phase  experiments  by 
Konrad  (1976,  s  -  1,  7)  are  also  shown.  Konrad  attempted  to  measure  the 
concentration  field  directly  in  dilution  experiments  in  non-homogeneous 
shear  layers.  While  he  used  the  term  entrainment  ratio,  his  definition 
involving  intermittency  was  an  attempt  to  included  the  mixed  fluid 
fraction  only,  and  as  a  result  will  be  called  here  the  mixed  fluid 
composition.  While  it  can  be  seen  that  this  assignement  of  a  value  to 
the  entrainment  ratio  follows  the  (geometrically  estimated)  Ev  trend,  it 
is  consistently  lower.  In  fact,  for  s  <  4,  the  mixed  fluid  ratio  EvM 
appears  to  be  a  nearly  fixed  fraction  of  the  geometric  entrainment  ratio 
^y,  i.e.  Ey^j/Ey  31  0.7  . 

There  is,  of  course,  no  a  priori  reason  to  assume  that  the  two 
quantities  should  be  the  same;  the  fluid  need  not  necessarily  be  mixed 
in  the  ratio  at  which  it  crosses  the  shear  layer  boundaries. 
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Nevertheless,  the  difference  in  the  behavior  for  density  ratios  greater 
than  unity  (low  speed  side  denser)  as  compared  to  density  ratios  smaller 
than  unity  (high  speed  side  denser)  is  noteworthy.  As  of  this  writing 
it  is  not  clear  whether  this  is  a  manifestation  of  the  different  kinds 
of  instability  modes  that  can  be  supported  by  shear  layers  at  different 
density  ratios  or  is  an  artifact  of  the  methods  we  are  using  to  estimate 
the  entrainment  ratio.  If  the  former  should  prove  to  be  the  case,  it 
would  constitute  a  new  result  in  the  fluid  mechanics  of  turbulent  shear 
layer  mixing. 

This  work  is  part  of  the  on-going  Ph.D.  thesis  work  of  Cliff 
Frieler . 


3.0  LIQUID  PHASE  SHEAR  LAYER 

This  work  has  been  completed.  An  early  announcement  in  the  form  of 
an  AIAA  Paper  was  published  during  the  present  reporting  period 
(Koochefahani  &  Diraotakis  1985).  A  more  extended  paper  documenting  the 
results  including  those  involving  chemical  reactions  has  been  accepted 
for  publication  in  the  of  Fluid  Mech.  (Koochesfahani  &  Dimotakis 
1986). 


4.0  TURBULENT  JET  MIXING  &  CHEMICAL  REACTIONS 


4 . 1  Gas  Phase  Jet  Mixing  -  Laser  Rayleigh  Scattering 


The  characterization  of  Schmidt  number  (ratio  of  kinematic 
viscosity  to  species  diffusivity)  effects  is  important  in  understanding 
turbulent  mixing  and  combusting  flows.  The  gas  phase  mixing  experiments 
described  here  are  directed  towards  determining  quantitatively  the 
differences  between  gases  and  liquids  for  turbulent  jets  and  plumes. 
Since  we  expect  that  the  effects  of  varying  diffusivity  would  be 
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confined  to  the  smallest  scales  of  the  flow,  high  resolution 
measurements,  in  both  time  and  space,  have  been  undertaken.  The 
development  of  a  laser  Rayleigh  scattering  diagnostic  for  these 
measurements  was  reviewed  in  the  interim  report  of  last  year.  This  year 
all  of  our  work  has  focused  on  the  full  scale  experiment  in  which  we 
will  make  the  high  resolution  measurements.  A  schematic  of  this 
facility  has  been  included  as  figure  5. 

The  construction  of  the  full  scale  experiment  facility  has  been 
completed.  This  includes  all  of  the  sub-systems  for  manipulation  of  the 
gases  and  optics.  Testing  of  these  systems  has  been  completed  with  good 
results  in  all  but  one  case.  A  few  preliminary  runs  were  made  as  a 
check  of  system  integration.  What  follows  in  this  section  is  a  summary 
of  the  completed  work. 

The  main  enclosure  on  which  the  laser  is  mounted  and  in  which  the 
gases  are  mixed  was  fabricated,  varnished,  painted,  sealed  and  placed  in 
a  laboratory  which  has  the  proper  electrical  capability  for  the  laser 
system.  Unfortunately,  this  was  in  the  sub-basement  and  it  was 
necessary  to  install  a  special  explosion  proof  suction  line  to  the  roof 
(five  stories  up)  in  order  to  safely  vent  the  exhaust  gases  from  the 
experiment.  The  large  optical  traverse  and  the  laser  beam  manipulation 
optics  were  installed  and  tested.  The  resolution  of  the  collection 
optics  was  extensively  tested  and  the  results  were  in  agreement  with  the 
calculated  values.  The  nozzle  assembly  was  completed  and  installed. 
The  exit  turbulence  level  was  found  to  be  less  than  0.2)1.  The  exit 
velocity  profile  was  found  to  be  uniform  to  within  the  accuracy  of  the 
pressure  transducer  (0.2$).  Combustible  gas  detectors  were  installed  in 
the  laboratory  as  a  safety  precaution  since  some  of  the  gases  we  plan  to 
use  are  volatile  (explosive)  hydrocarbons .  Appropriate  data  acquisition 
and  processing  software  was  developed.  The  construction  of  the  gas 
delivery  systems  has  been  completed.  Their  testing  has  not  been 
completed  since  the  required  long  time  stability  of  the  flows  from  these 
systems  was  not  achieved.  The  gas  delivery  system  for  the  jet  was 
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modified  to  include  a  driver  vessel  which  would  assure  a  constant 
driving  pressure  which  was  unattainable  by  an  earlier  design.  The 
temperature  of  the  gas  delivered  by  the  coflow  system  was  found  to  drop 
by  about  1.5  K  during  one  minute.  Unfortunately  this  temperature 
difference,  in  addition  to  the  density  differences  resulting  from  air 
left  in  the  main  enclosure  from  inadequate  purging,  was  sufficient  to 
cause  unsteady  buoyancy  currents  in  the  main  enclosure  after  about  15 
seconds  of  run  when  argon  was  used  as  the  reservoir  gas.  This  problem 
is  being  addressed  with  a  feedback  heating  system  and  a  new  purging 
procedure. 

Since  the  troubles  with  the  gas  delivery  system  are  only  a 
consideration  affecting  longer  run  times,  some  short  preliminary  runs 
were  made.  This  was  a  useful  test  of  all  aspects  of  the  entire 
experiment  and  has  yielded  some  data  that  already  suggest  qualitative 
differences  between  gas  phase  and  liquid  phase  (Sc  »  1 )  mixing.  In 
particular,  the  gas  phase  data  suggest  that  reservoir  fluid  is  not  found 
on  the  centerline  of  the  jet  as  frequently  as  in  liquids  and  gradients 
in  the  concentration  field  of  the  gas  phase  jet  are  seen  to  be  smaller 
for  the  same  Reynolds  number.  We  should  repeat  here  that  the  temporal 
and  spatial  resolution  in  these  measurements  is  adequate  to  resolve  the 
theoretically  computed  corresponding  scales  in  the  flow.  A  time  trace 
of  some  data  from  these  runs  has  been  included  a3  figure  6. 

The  work  planned  for  next  year  is  to  correct  the  remaining  minor 
problems  with  the  gas  delivery  systems  and  make  the  detailed 
measurements  that  will  yield  quantitative  information  about  the 
differences  between  gas  and  liquid  phase  turbulent  jet  mixing. 

This  work  is  part  of  the  on-going  Ph.D.  research  by  David  Dowling. 
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4.2  Liquid  Phase  Jet  Mixing  &  Chemical  Reactions 


The  liquid  phase  jet  mixing  experiments  were  completed  during  this 
reporting  period.  The  results  have  been  documented  in  the  form  of  a 
GALCIT  Ph.D.  thesis  (Dahm  1985).  A  first  paper  of  aspects  of  these 
data  that  was  presented  earlier  (Dahm  &  Dimotakis  1985)  is  in  the  final 
stages  of  preparation  to  be  submitted  for  publication  in  the  AIAA 
Journal . 


4.3  Soot  scattering  measurements  in  turbulent  flames 

We  are  finishing  a  report  (Miake-Lye  &  Toner  1986)  on  the 
experimental  soot-scattering  investigation  carried  out  during  the  the 
past  year  on  a  large  turbulent  buoyant  diffusion  flame.  The 
experimental  apparatus,  the  soot-scattering  data,  and  the  data  analysis 
procedures  are  described  in  in  this  report.  Several  important 
conclusions  can  be  drawn  about  the  soot  distribution  in  the  large 
buoyant  diffusion  (methane)  flame  we  have  been  studying  at  the  initial 
phases  of  this  effort. 

In  our  report,  we  have  set  forth  arguments  justifying  the  use  of 
soot  as  a  marker  for  the  combustion  interface  Ln  the  flame  under  study . 
These  arguments  are  based  on  the  facts  that 

1 .  essentially  all  the  soot  generated  by  the  flame  is  consumed  by 
the  flame  tip  and,  further, 

2.  there  is  no  evidence  for  cold  soot  that  is  formed  at  one  stage 
and  consumed  at  a  later  stage. 

These  data,  in  conjunction  with  our  r-t  images,  lead  us  to  conclude  that 
the  scot  is  localized  in  a  thin  layer  close  to  and  locally  following  the 
reaction  interface  surface. 
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Thus,  our  imaging  measurements  of  the  soot  distribution,  wherein  we 
have  imaged  a  line  segment  within  the  flame  as  a  function  of  time,  allow 
us  to  make  inferences  about  the  nature  of  the  topology  of  the  combustion 
zone  itself.  Construction  of  r-t  images  from  many  individual  scans  of 
the  imaged  line  segment  indicates  that  the  soot  is  found  in  relatively 
thin,  contorted  zones  which  we  interpret  to  be  convoluted  laminar 
diffusion  flamelets.  The  flamelets  are  continuous  near  the  burner,  but 
are  broken  and  discontinuous  further  downstream  in  the  flame  where  the 
velocities  are  greater. 

Based  on  our  r-t  data,  we  have  calculated  the  soot  intermittency , 
or  fraction  of  time  that  soot  is  present,  at  each  radial  station.  The 
soot  intermittency  is  highest  near  the  center  and  falls  to  zero  at  the 
edges  of  the  flame.  By  determining  the  instantaneous  flame  boundary,  we 
have  also  calculated  the  conditional  soot  intermittency,  which  we  define 
to  be  the  fraction  of  time  that  soot  is  present  within  the  flame 
boundary  at  any  radial  station.  The  conditional  soot  intermittency  is 
nearly  constant  across  the  flame.  These  conclusions  hold  true  for  the 
various  axial  stations  and  the  range  of  fuel  flow  rates  that  we  have 
investigated. 

Our  measurements  to  date  lead  us  to  a  picture  of  the  flame  made  up 
of  thin  laminar  diffusion  flamelets  distributed  fairly  uniformly  within 
the  large  puff-like  envelope  of  the  flame.  The  interior  of  the  puff  is 
entirely  free  of  soot  for  40  to  60?  of  the  time;  the  remainder  of  the 
time  some  part  of  a  relatively  thin  soot-sheet  is  present.  Since  the 
width  of  the  reaction  zone  is  probably  smaller  than  the  width  of  the 
soot-sheet  and  since  soot  may  strain  out  before  the  reaction  is 
quenched,  the  fraction  of  time  that  the  reaction  interface  is  present 
may  be  different. 


Investigations  to  follow  will  make  use  of  a  new  burner  designed  to 
permit  the  illuminating  laser  beam  to  be  directed  parallel  to  the  burner 
axis  at  several  radial  stations.  The  burner  construction  is  now 
complete  and  we  are  integrating  it  into  our  experimental  apparatus.  By 
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permitting  tthe  laser  beam  to  be  directed  parallel  to  the  axis  of  the 
burner,  this  burner  will  allow  us  to  measure  the  celerity  of  the  puff 
structures.  Modifications  of  our  optical  geometry  are  presently  in 
progress . 

We  will  also  make  higher  resolution  1 — t  measurements  beyond  x/D  of 
0.75  to  complement  our  work  to  date.  We  will  have  to  observe  a  smaller 
section  of  the  flame  in  order  to  resolve  the  smallest  soot-sheets  but, 
since  our  previous  data  include  images  of  the  entire  flame  diameter,  we 
will  be  able  to  piece  together  a  more  complete  view  of  the  soot 
structure  even  at  these  axial  stations.  These  experiments  will  make  use 
of  the  existing  experimental  apparatus  and  recent  improvements  in  our 
real-time  digital  image  acquisition  capabilities.  Subsequent,  planned 
for  experiments  will  push  the  flow  from  the  buoyancy-driven  regime  into 
one  that  is  more  closely  momentum  driven. 

These  studies  are  co-sponsored  by  the  Gas  Research  Institute  (GRI). 

5.0  COMPUTATIONAL  EFFORT 

The  objective  of  this  research  is  to  numerically  simulate  a  2-D 
incompressible  shear  layer  with  the  capability  of  including  chemical 
reactions.  By  Kelvin  and  Helmholtz's  theorems,  the  inviscid  motion  of 
vorticity  at  high  Reynolds  numbers  is  determined  by  the  local  fluid 
velocity  which  is  itself  based  on  the  vorticity  field.  The  techniques 
of  discretizing  vorticity  and  calculating  the  mutual  interactions  have 
been  frequently  used  to  simulate  the  flow.  The  point  vortex  and  vortex 
blob  are  two  typical  methods  of  discretization. 


k**-**-**-^-  '*a~ J'-  *ALj&iU  *-  -V*v-  V-  M-  A  V-  jU  »•-  W-  n.  V-  amsC*.  1 kVjtV  a\ 


The  effort  so  far  has  mainly  concentrated  on  improving  the 
computational  efficiency  and  accuracy.  There  are  several  ways  to 
accomplish  this  purpose.  One  of  the  methods  of  improving  the 
computational  efficiency  is  to  use  special  algorithms  which  can  reduce 
the  calculations  of  the  vortex  interactions  from  O(N^)  to  O(NlogN), 
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where  N  is  the  number  of  computational  elements  (vortices).  An 
independent  method  (which  can  be  used  in  conjunction  with  the  N2  to 
N  logN  reduction)  is  to  modify  the  blob  method  itself. 

We  have  been  investigating  the  possibility  of  using  an  elliptical 
guassian  vortex  blob  instead  of  the  regular  circular  guassian  vortex 
blob,  in  the  hope  that  the  number  of  vortices,  required  for  the  same 
accuracy,  can  be  reduced  to  model  the  same  flow.  The  reason  for  this  is 
as  follows.  Since  a  concentration  of  vorticity  causes  a  strain  field 
around  itself,  the  effects  on  the  motion  of  the  other  vortices  from  that 
strain  field  should  be  considered.  The  analytical  result  by  Kida  (1981) 
shows  that  an  elliptic  vortex  of  uniform  vorticity  in  a  uniform  strain 
field  rotates  around  its  center,  preserves  the  elliptical  shape  and 
keeps  the  area  constant.  However,  the  rates  of  changes  in  the  aspect 
ratio  and  rotation  angle  are  not  constant.  See  figure  7.  Consequently, 
to  model  the  fluid  precisely,  the  blobs  should  be  allowed  to  deform  in 
response  to  the  strain  fields  from  the  other  vortices. 

By  analysis,  the  motion  of  the  elliptical  guassian  blob  vortex  in 
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translation  with  the 
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by 

vortex 
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interaction 

y. 

£ 

n 

2. 

3elf-induced  rotation 

R 

3. 

deformation  caused  by 

strain 

fields. 

It  can  also  be  shown  that  the  motion  of  each  blob  depends  only  on 
the  present  value  of  the  strain  fields  and  the  state  of  the  blob  itself. 
The  rate  of  deformation  can  be  obtained  analytically  if  the  velocity 
field  around  the  blobs  are  known  and  those  effects  from  the  strain 
fields  are  linearly  summed  up.  See  figure  8. 
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We  are  currently  studying  the  self-induced  rotation  and  the 
velocity  field  around  an  elliptical  guassian  core.  Unfortunately,  exact 
solutions  are  not  available.  However,  the  self-induced  rotation  of  the 
blob  is  asymptotic  to  that  of  the  uniform  vorticity  and  the  velocity 
field  may  be  obtained  by  multipole  expansion  methods.  It  is  clear  that 
this  new  blob  method  needs  more  calculations  and  bookkeeping  per  vortex 
blob  than  the  regular  blob  method.  It  is  anticipated,  however,  that  the 
total  CPU  time  will  be  less  as  a  result  of  the  reduction,  for  a  given 
required  accuracy,  of  the  required  number  of  vortices. 

This  effort  is  part  of  the  on-going  Ph.D.  thesis  work  of  Kayo  Ide. 
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ABSTRACT 

The  linear  spatial  instability  of  two-dimensional  two-stream  plane 
mixing  layers  has  been  studied  extensively  in  the  past.  In  the  case  of 
uniform  density,  Michalke  (1965)  investigated  the  single-stream  shear 
layer  and  Monkewitz  &  Huerre  (1982)  considered  the  effect  of  the 
velocity  ratio.  Maslowe  &  Kelly  (1971)  studied  the  stratified 
(non-uniform  density)  shear  layers  and  snowed  that  density  variations 
can  be  destabilizing.  In  all  these  studies,  the  mean  velocity  profile 
has  been  assumed  to  be  monotonically  increasing  from  the  value  on  the 
low-speed  stream  to  that  on  the  high-speed  stream  and  usually  the 
hyperbolic  tangent  form  is  used.  It  should  be  noted,  however,  that 
under  experimental  conditions  the  initial  mean  velocity  profile  almost 
always  has  a  wake  component  due  to  the  boundary  layers  on  the  two  sides 
of  the  splitter  plate.  The  effect  of  the  wake  component  has  only 
recently  come  into  consideration  with  the  investigations  of  Miau  1984 
and  Zhang  et  al.  1984  for  the  uniform  density  case. 

The  purpose  of  the  present  work  is  to  study  the  instability 
characteristics  of  both  uniform  and  non-uniform  density  plane  shear 
layers  taking  into  account  the  wake  component  of  the  initial  velocity 
profile.  The  inviscid,  linear,  parallel-flow  stability  analysis  of 
spatially  growing  disturbances  is  utilized  to  numerically  calculate  the 
range  of  unstable  frequencies  and  wave-numbers. 

We  consider  the  general  case  of  a  two-stream  plane  shear  layer  with 
Uj,  pt  as  the  free-stream  velocity  and  density  on  the  high-speed  stream 
and  U2 ,  p2  as  the  corresponding  quantities  on  the  low-speed  side  of  the 
layer.  All  quantities  used  here  are  normalized  with  the  average 
velocity  (Ul+U2)/2,  average  density  (p,+p2)/2  and  the  local  layer 
thickness  <5  as  the  length  scale.  The  mean  velocity  profile  is  composed 
of  the  usual  hyperbolic  tangent  profile  plus  a  wake  component 
represented  by  a  Gaussian  distribution  and  has  the  form 

U(7)  -  1  "  Au  tanh(Y)  -  W  e'y2 


2 


where  Au  »  (U  l -U2 )  /  ( U ,  +U2 )  and  W  is  the  normalized  wake  deficit.  The 
mean  density  profile  has  a  hyperbolic  tangent  profile  and  is  given  by 

p(Y)  -  1  +  Ap  tanh[ ( Y-Y 0 )/o] 

where  Ap  -  ( p l -p2 ) / ( p t +p2 )  and  Ya  and  o  adjust  the  lateral  position  and 
thickness  of  the  density  profile  relative  to  the  velocity  profile. 

The  disturbance  stream  function  is  written  in  the  form 

^  =  $ ( Y)  e  i ( -  Bt ) 

where  a  -  ar+ia^  is  the  complex  wave-number  and  8  is  frequency  which  is 
real  for  the  spatial  case.  For  the  case  of  incompressible  flow  and 
negligible  buoyancy  effects  (i.e.  gravity  is  ignored),  it  can  be  shown 
that  the  disturbance  eigenfunction  $  satisfies  the  equation 


,  „  U"  +  p’  U’/p 

t>"  +  (p'/p)  <f> '  -  [a2  +  -  1 

U-B/a 


where  (  )'  corresponds  to  d/dY.  The  equation  above  reduces  to  the 
Rayleigh  equation  when  the  density  is  uniform.  A  "shooting"  technique 
is  used  to  solve  this  eigenvalue  equation  and  calculate  the  spatial 
growth  rate,  -a^ ,  of  unstable  disturbances  and  the  corresponding 
wave-number,  ar ,  versus  frequency  3. 


Representative  results  for  the  uniform  density  case  are  shown  in 
figure  1  .  The  main  result  is  that  when  the  wake  component  is  present 
there  are  two  unstable  branches  as  opposed  to  one  in  the  case  of  the 
hyperbolic  tangent.  These  two  branches  do  not  have  much  in  common  with 
the  "sinuous"  and  "varicose"  modes  of  instability  calculated  by 
Mattingly  &  Criminale  (1972)  in  the  case  of  pure  wakes  since  the 
symmetry  conditions  that  exist  in  the  pure  wake  do  not  hold  here.  The 
stronger  branch  leads  to  the  usual  Kelvin-Helmholtz  roll-up  patterns 
observed  in  shear  layers  and  it  is  referred  to  here  as  the  shear  layer 
branch.  The  lower  branch  is  suspected  to  correspond  to  roll-up  patterns 


that  resemble  those  in  pure  wake  flows  and,  therefore,  we  call  it  the 
wake  branch.  This  branch,  in  the  uniform  density  case,  is  rarely 
observed  since  its  amplification  rate  is  quite  lower  than  the  shear 
layer  branch. 

To  investigate  the  effect  of  non-uniform  density,  we  consider  the 
case  of  low-speed  stream  having  the  higher  density  (see  figure  2).  The 
density  profile  has  been  arranged  to  have  its  inflection  point  at  the 
minimum  of  the  velocity  profile  and  its  thickness  much  thinner  than  that 
of  the  velocity  profile.  These  conditions  are  expected  to  hold  in  the 
initial  region  of  the  flow  near  the  splitter  plate.  The  qualitative 
features  of  the  results  are  not  sensitive  to  these  conditions  as  long  as 
the  density  profile  is  reasonably  thin  relative  to  the  velocity  profile. 

The  results  for  the  case  of  non-uniform  density  are  shown  in 
figure  3.  In  figure  3a,  the  uniform  density  case,  the  two  instability 
branches  similar  to  figure  1b  can  be  seen.  The  most  important  result  is 
that  when  the  high  density  is  on  the  low-speed  side,  the  two  branches 
have  similar  growth  rates.  In  fact,  the  normally  weak  branch  (wake 
branch)  in  the  uniform  density  case,  now  seems  to  have  a  slightly  larger 
growth  rate  than  the  shear  layer  branch.  This  means  that,  depending  on 
the  spectrum  of  disturbances  in  the  flow  and  the  extent  of  the 
persistence  of  the  wake  component  in  the  downstream  region,  a  shear 
layer  of  non-uniform  density  may  not  roll-up  like  the  usual 
Kelvin-Helmholtz  structures  but  more  like  a  wake.  The  finding  in 
figure  3  also  suggests  that,  under  the  right  flow  conditions,  both  shear 
layer  and  wake  modes  of  instability  may  exist  simultaneously  and 
interact  with  each  other. 

To  show  that  the  features  described  above  are,  in  fact,  possible, 
the  shear  layer  between  a  high-speed  stream  of  Helium  and  low-speed 
stream  of  Argon  was  forced  acoustically.  The  flow  visualization  by 
Schlieren  photography,  figure  4,  shows  that  both  wake  and  shear  layer 
modes  can  be  generated  in  a  two-stream  shear  layer. 
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Figure  3a.  Instability  characteristics 
of  the  mean  velocity  profile 
of  figure  2  with  uniform 
density. 
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Figure  3b.  Instability  characteristics 
of  the  mean  velocity  and 
density  profiles  of  figure 
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Figure  4 


Schlieren  photographs  of  the  shear  layer  mode  (a) 
and  the  wake  mode  (b).  Flow  is  from  left  to  right 
with  high-speed  stream  on  top  (l^/U^  =*  0.38, 
P2/P1=*10). 
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ABSTRACT 

A  chemical  reaction  for  which  the  reaction  rate  can  be  varied  Is  studied 
in  a  fully  developed,  two-dimensional,  turbulent  mixing  layer.  The  layer 
is  formed  between  two  nitrogen  streams,  one  carrying  low  concentrations  of 
fluorine  and  the  other  hydrogen  and  nitric  oxide.  For  fixed  concentrations 
of  fluorine  and  hydrogen  and  for  nitric  oxide  concentrations  that  are 
small  fractions  of  the  fluorine  concentration,  the  heat  release  is  fixed 
but  the  overall  reaction  rate  is  controlled  by  the  nitric  oxide  concentra¬ 
tion.  Therefore,  for  fixed  flow  conditions,  the  nitric  oxide  concentra¬ 
tion  determines  the  ratio  of  the  reaction  rate  to  the  mixing  rate.  For 
large  values  of  this  ratio,  the  amount  of  product,  at  a  given  downstream 
location,  measured  by  the  mean  temperature  rise,  is  independent  of  the 
reaction  rate,  i.e.,  the  reaction  is  mixing  limited.  As  the  reaction  rate 
is  reduced  the  major  effects  are:  (1)  amount  of  product  declines  (as 
expected),  (2)  the  mean  temperature  profile,  which  is  initially  somewhat 
unsymmetrlcal  because  the  hydrogen-fluorine  freestream  concentration  ratio 
is  set  at  a  large  value,  becomes  symmetrical,  and  (3)  the  ramp-like 
instantaneous  temperature  traces  within  the  large  structure  gradually 
become  more  top-hat. 
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1 .  Introduction 

This  work  represents  the  third  in  a  series  of  studies  aimed  at 
achieving  an  understanding  of  mixing  and  combustion  in  a  fully  developed, 
two-dimensional,  turbulent  mixing  layer,  Fig.  1.  In  these  experiments, 
the  reactants  are  low  concentrations  of  hydrogen  and  fluorine  carried  in 
separate  nitrogen  freestreams.  Nitric  oxide,  in  concentrations  that  are 
small  compared  to  the  fluorine,  is  carried  with  the  hydrogen  stream.  When 
the  two  streams  mix,  the  nitric  oxide  dissociates  a  fraction  of  the 
fluorine  to  initiate  and  control  the  hydrogen-fluorine  chain  reaction.  In 
all  experiments,  the  heat  release  is  sufficiently  low  that  the  heat 
release  does  not  couple  with  the  fluid  mechanics,  so  that  the  overall  pro¬ 
perties  of  the  mixing  layer  are  not  significantly  changed  from  those  of 
the  non-reacting  case.  It  is  important  to  note  that  the  present  chemical 
system  is  hypergolic  (low  activation  energy)  at  all  concentrations,  unlike 
hydrocarbon  -  air  systems  which  require  ignition  sources  to  sustain 
combustion,  have  flammability  limits  and  can  be  quenched  by  high  values  of 
strain  rate. 

The  first  of  these  efforts  [1,2]  examined  the  effects  of  the 
hydrogen-fluorine  freestream  concentration  ratio  upon  the  amount  of  pro¬ 
duct  formed  in  the  fully  developed  layer  at  a  single  Reynolds  number  of 
4*  10®  (based  on  velocity  difference  and  downstream  distance).  The 

results  confirmed  the  central  role  of  the  large  structure  dynamics  in  mix¬ 
ing  and  combustion.  The  asymmetry  in  entrainment  of  the  mixing  layer  was 
demonstrated  through  'flip'  experiments  whereby  the  lean  reactant  being 


carried  on  the  high-speed  or  low-speed  side  of  the  layer  leads  to  dif¬ 
ferent  amounts  of  product.  The  experiments  also  showed  that  when  one 
reactant  concentration  is  fixed,  a  nearly  asymptotic  product  profile  is 
reached  when  the  other  reactant  concentration  is  approximately  eight  times 
the  first.  Under  these  conditions,  essentially  all  the  molecularly  mixed 
lean  reactant  is  completely  consumed.  It  is  instructive  to  note  that  for 
chemical  reactions  in  water  at  similar  Reynolds  numbers,  Koochesf ahani  [3] 
found  that  there  was  approximately  50%  less  (normalized)  product  than  in 
the  gas.  This  has  led  to  the  conclusion  that  the  Schmidt  number  influ¬ 
ences  the  amount  of  product  formed  in  turbulent  flows,  even  at  Reynolds 
numbers  as  high  as  4*10®. 

In  the  second  effort  [4],  the  dependence  of  the  amount  of  product 
formed  in  the  mixing  layer  upon  Reynolds  number  was  studied.  It  was  found 
that  the  large  organized  structures  are  present  at  all  Reynolds  numbers 
investigated,  irrespective  of  whether  the  boundary  layer  on  the  splitter 
plate  is  laminar  or  turbulent,  and  that  the  product  dependence  on  Reynolds 
number  is  weak.  There  is,  approximately,  a  20%  decrease  in  the  amount  of 
product  formed  for  a  factor  of  10  increase  in  Reynolds  number  (or  6% 
decrease  per  factor  of  2  in  Reynolds  number). 

In  both  of  the  above  investigations,  the  amount  of  nitric  oxide  was 
sufficiently  large  to  insure  that  the  reaction  was,  to  a  good  approxima¬ 
tion,  mixing  limited.  In  the  present  investigation  we  reduce  the  nitric 
oxide  concentration  to  study  the  effects  of  the  chemical  reaction  rate.  A 
chemical  kinetics  code,  CHEMKIN  [5],  was  used  to  study  the  hydrogen- 
fluorine-nitric  oxide  system.  The  resuits  show  that  the  nitric  oxide  con¬ 
centration  controls  the  overall  reaction  rate,  while  the  heat  of  reaction 
is  the  same.  Since  the  runs  described  below  are  all  performed  at  low  heat, 
the  fluid  mechanics  also  remains  unchanged. 

In  the  investigations  mentioned  above,  the  findings  were  interpreted 
in  the  context  of  the  simple  model  for  mix.'ng  and  product  formation  in  the 
turbulent  shear  layer  proposed  by  Broadwell  4  Breidenthal  [6j.  The 


present  results  are  interpreted  in  terms  of  this  model  in  Broadwell  4  Mun- 
gal  [7], 


2.  The  Chemical  Reactions 

The  principal  reactions  in  the  experiment  consists  of  the  hydrogen- 
fluorine  chain  reaction  and  the  nitric  oxide  fluorine  dissociation  reac¬ 
tion  . 
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where  k  is  given  in  cc/mole-s,  T  in  degrees  K  and  R  is  the  universal  gas 
constant  in  cal/mole-K. 


Reactions  (1)  and  (2)  are  the  so-called  cold  and  hot  reactions, 
respectively.  The  cold  reaction,  which  carries  25%  of  the  total  heat 
release,  is  faster  than  the  hot  reaction  by  about  an  order  of  magnitude  at 
300K  (kj^  :  i.exio13,  k2  :  9.5X1011)-  The  rate  constants  ^  and  k2  are 
taken  from  Cohen  4  Bott  [8]  and  k3  from  Baulch  et  al .  [9], 

Additional  reactions  which  must  be  considered  are: 
k4 

F  -*•  NO  +  M  ♦  NOF  +  M,  AQ  =  -57  kcal  #  ^  ;  3xl016  (4) 


and 


m 


Fjtmm 

mMJ 


ISS 


r 


-  4  - 


F 

F 

*  M  - 

F2 

*  M, 

k5  : 

14 

3.2x10 

(5) 

k6 

2.5*1016 

H 

*  F 

+  M  - 

HF 

+  M, 

k6  : 

(6) 

k? 

3.0xl015 

H 

*  H 

M  -* 

H2 

+  M, 

k7  " 

(7) 

where  the  k4  _  rates  are  given  in  cc2/»ole2-s  and  are  estimated  at  300K 
from  Baulch  et  al.  [9]. 

We  have  used  the  CHEMKIN  [5]  computer  code  to  study  the  reaction  set 
(1)  -  (7)  in  a  premixed  constant  pressure  reactor.  The  initial  condition 

’  0.5%Fg,  0.015%N0  in  a  Ng  diluent  simulate  our  normal  operating  con¬ 
ditions  of  S'Hg,  0.03%NO  in  Ng  on  the  high-speed  side  and  l%Fg  in  Ng  on 
the  low-speed  side  (this  nominal  condition  will  be  referred  to  as  [NO]*). 
The  CHEMKIN  code  computes  the  species  concentration  and  temperature  as  a 
function  of  time  for  the  mixture  at  constant  pressure  or  constant  volume, 
the  former  being  the  option  applicable  to  this  work.  The  result  of  such  a 
calculation  is  shown  in  Fig.  2.  The  starting  temperature  is  300K  and  the 
system  asymptotes  to  a  temperature  of  393K.  The  figure  shows  the  results 
if  only  the  first  three  reactions  ( 1 >  —  ( 3 )  are  used,  together  with  that  for 
all  seven  reactions  (l)-(7).  The  agreement  suggests  that  reactions  ( 1 ) - 
(3)  determine  the  behaviour  of  the  chemical  system  at  these  temperatures, 
while  ( 4 ) - ( 7 )  are  relatively  less  important.  We  have  fitted  a  straight 
line  to  the  curve  shown  in  Fig.  2  and  taken  the  time  required  for  a  tem¬ 
perature  rise  from  300K  to  393K  to  be  a  measure  of  the  chemical  time  of 
the  system.  Fig.  2  shows  that  for  the  rate  constants  shown  above,  the 
chemical  time  for  our  nominal  condition,  t*.  is  0.53  ms.  (It  is 
worthwhile  to  note  that  if  we  simulate  an  operating  condition  of  1%K2 

0.03%  NO  together  with  1%F^  in  a  Ng  diluent,  as  was  performed  in  Reference 
2.  the  chemical  time  Increases  by  about  15%. ) 

Table  1  shows  (among  other  results)  the  overall  chemical  rate,  k  (the 
inverse  of  the  chemical  time,  t)  as  computed  by  CHEMKIN,  for  various  con¬ 
centrations  of  nitric  oxide  for  the  reaction  set  ( 1 )  - ( 7 )  .  As  is  shown 
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analytically  in  Appendix  1,  the  chemical  rate  is  nearly  proportional  to 
[NO]  at  very  low  concentrations  (relative  to  [F2])  and  varies  as  the 
square  root  of  [NO]  at  higher  concentrations.  In  summary,  the  chemical 
system  (l)-(7)  describes,  for  our  conditions,  a  reaction  in  which  the  heat 
release  is  fixed  by  the  hydrogen  and  fluorine  concentrations  but  the 
overall  rate  is  controlled  by  the  nitric  oxide  concentration. 

3.  Experimental  Apparatus 

The  experimental  facility  has  been  described  previously  [1,2]  and 
will  only  be  discussed  briefly  here.  It  is  of  the  blowdown  type  in  which 
premixed  volumes  of  hydrogen  and  nitric  oxide  in  nitrogen,  and  fluorine  in 
nitrogen  are  discharged  through  sonic  orifices,  maintaining  a  constant 
mass  flux  in  each  of  the  freestreams.  Both  streams  enter  a  settling  and 
contraction  section  for  turbulence  reduction  with  the  high  speed  flow 
emerging  from  a  6:1  contraction  at  an  exit  area  of  5  x  20  cm,  and  the  low 
speed  emerging  from  a  4:1  contraction  in  a  7.5  x  20  cm  exit  area  The 
flows  then  meet  at  the  tip  of  a  splitter  plate  as  shown  in  Fig.  1.  All 
measurements  are  recorded  at  a  station  x  =  45.7  cm  downstream  of  the 
splitter  plate  trailing  edge. 

Runs  were  performed  for  flow  velocities  of  =  22  and  44  ra/s  with  a 
fixed  speed  ratio  =  0.40.  The  Reynolds  numbers  based  on  the  velo¬ 
city  difference  AU,  and  the  1%  thickness  of  the  layer  <5^  (distance  between 
points  at  which  the  mean  temperature  rise  is  1%  of  the  maximum  mean  tem¬ 
perature  rise),  and  the  cold  kinematic  viscosity  of  nitrogen  are  respec¬ 
tively  6.7xl04,  and  1.3xl05.  These  values  are  all  beyond  the  mixing  tran¬ 
sition  observed  by  Konrad  [10]  and  Breidenthal  [11].  The  high  speed  tur¬ 
bulence  level  is  about  0 . 7\  rms  at  the  lower  speed  and  decreases  somewhat 
at  the  higher  3peed.  For  =.  22  m/s,  the  high-speed  boundary  layer 
momentum  thickness.  was  estimated  using  Thwaites'  method,  with  the 
Reynolds  number  v  estimated  to  be  240.  Thus  at  the  measuring  sta¬ 
tion,  x/®1  ;  2800. 


4.  Measurements 


For  each  run,  the  streamwise  pressure  gradient  is  set  to  zero  by 
adjustment  of  the  low  speed  sidewall.  Temperature  rise  is  recorded  by  a 
rake  of  8  cold  wire  resistance  thermometers  placed  across  the  width  of  the 
layer,  each  driven  by  a  constant  current  of  0.4mA.  The  wires  are  made  of 
2.5  Mm  diameter  platinum  /  10%  rhodium  welded  to  inconel  prongs  with  a 
span  of  1.5  mm.  The  ratio  of  the  length  of  the  wire  to  Betchov's  cold 
length  [12]  is  about  20  at  the  lower  Reynolds  number.  During  a  run,  each 
wire  is  sampled  at  10  and  12.5  kHz  for  the  two  speeds,  corresponding  to 
total  data  rates  of  80  and  100  kHz. 

k 

As  described  in  [2],  the  cold  wires  are  calibrated  directly  before 
each  run  and  record  an  accurate  mean  temperature  profile,  while  excursions 
from  the  mean  suffer  from  conduction  error  due  to  the  prongs  and  the  ther¬ 
mal  lag  of  the  wire  itself  [13],  [14].  The  mean  temperature  profiles  were 
spot  checked  by  a  0.25  mm  diameter  chromel-alumel  thermocouple  at  a  single 
point  on  the  mean  temperature  rise  profile,  with  the  agreement  being  gen¬ 
erally  within  3%  or  better.  The  mean  temperature  profiles  shown  below  are 
repeatable  from  run  to  run  to  within  2  percent.  The  absolute  reactant 
concentrations  however,  are  known  only  to  approximately  3-5%,  and  thus  the 
absolute  value  of  the  potential  heat  release  has  an  uncertainty  of  the 
same  order. 


5.  Experimental  Results  and  Discussion 

For  the  runs  reported  here,  the  high-speed  stream  consisted  of  8%H2 
92%N2  together  with  varying  small  concentrations  of  NO,  while  the  low- 
speed  stream  consisted  of  0.96%F2  and  99%N2  (passivation  of  upstream 
screens  and  plates  removed  an  estimated  0.04%F2  froB  the  nominal  value  of 
I'Fg)-  with  an  adiabatic  flame  temperature  rise  of  159K.  As  was  noted 
above,  this  choice  of  reactant  concentrations  guarantees  that,  when  the 
reaction  rate  is  sufficiently  high,  the  fluorine  (lean  reactant)  which  is 


molecularly  mixed  with  the  hydrogen  is  aimost  totally  consumed. 


Figure  3a  shows  the  mean  temperature  rise  profiles  obtained  at  the 

x  ~  xq  =  45.7  cm  station  for  the  cases  of  decreasing  nitric  oxide  concen¬ 
tration.  namely  [NO]/[NO]*  =  l,  1/2,  1/4,  1/8.  1/16,  1/23.  1/32.  0  at 

=  22  m/s  (note  that  the  accuracy  in  loading  our  nominal  concentration 
of  NO  is  probably  within  3%  but  becomes  progressively  worse  as  the  amount 
of  NO  decreases).  The  symbols  indicate  the  measured  average  temperature 
rise  at  the  positions  of  the  eight  probes,  and  represents  the  average  of 
12,288  points  per  probe  during  which  time  at  least  100  structures  would 
have  passed  the  measuring  station.  The  smooth  curve  is  an  exponential  fit 
of  the  type 


T(ij)/Tf lm  =  exp  (  cQ  +  Cjj  *  Cgj2  +  c^3  +  c4j4  ) 

where  T  is  the  local  mean  temperature  rise,  Tflm  is  the  adiabatic  flame 

temperature  rise.  Jf  a  y/(x-Xg)(  ancj  the  point  v  *  0  corresponds  to  the 
position  of  the  dividing  streamline  (see  Table  3). 

We  see  that,  as  [NO]  decreases,  there  is  a  gradual  loss  in  area  under 
each  of  the  product  profiles,  together  with  a  tendency  for  the  profile  to 
become  more  symmetric.  In  fact,  for  the  case  [N0]/[N0]*  =  1/32,  the  pro¬ 
file  is  remarkably  symmetric  in  spite  of  the  fact  that  the  high-speed 
reactant  is  eight  times  more  concentrated  than  the  low-speed  reactant. 
This  trend  towards  symmetry  is  shown  more  directly  in  Fig.  3b  where  ail 
profiles  have  been  scaled  to  the  same  height  and  Fig.  3c  where  the  cases 

[ NO ] / [ NO ]  =  l,  1/32  are  shown.  It  should  also  be  noted  that  the  case 

* 

[ NO ] / [ NO ]  =  0  shows  no  product  both  experimentally  and  in  computation 

with  the  CHEMKIN  code. 

Figure  4  shows  the  mean  temperature  rise  profiles  for  the  cases  of 
increasing  the  nitric  oxide  concentration,  namely  [N0l/[N0]*  -  1,  3/2.  2. 
Table  1  indicates  that  the  effective  reaction  rates  for  the  latter  two 
runs  are  larger  than  the  nominal  by  factors  of  1.3  and  1.5  respectively. 
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The  similarity  of  the  profiles  implies  that  the  reaction  is  mixing  rate 
limited  for  [NO)  at  the  nominal  value  and  above. 


Figure  5a  shows  similar  results  taken  at  the  same  measuring  station, 
but  with  =  44  m/3  (and  Ug  =  0.4U^  to  maintain  the  same  speed  ratio)  for 
the  cases  [NO]/[NO]  »  2,  1,  1/4,  1/16  (see  also  Table  2.4).  The  scaled 
profiles  for  the  cases  [NO]/[NOj*  =  1.  1/16  are  shown  in  Fig.  5b  to  again 
demonstrate  the  symmetrizing  of  the  mean  temperature  profile.  It  is 
worthwhile  to  note  that,  in  the  fast  chemistry  limit,  the  absolute  area 
under  the  temperature  profile  for  the  =  44  „/s  case  iS  less,  by  about 
10  percent,  than  that  for  Uj  =  22  m/s.  As  is  discussed  in  [4],  this  is 
believed  to  be  a  Reynolds  number  effect  and  not  a  Oamkohler  number  effect, 
since  increasing  the  [NO]  concentration  by  a  factor  of  2  at  both  Reynolds 
numbers  with  the  corresponding  increase  in  chemical  rates,  does  not  elim¬ 
inate  the  difference. 

The  amount  of  product  generated  by  the  reaction  is  characterized  by 
the  product  thickness 


4 

P 


«  cp(y) 

*  /  -I— 


00 

/ 

—00 


cD  t(v> 

coo  AQ 


dy 


where  cp(y)  is  the  concentration  of  product  and  c,,,  is  the  freestream 
fluorine  concentration.  Since  product  is  analogous  to  temperature  rise, 
the  second  expression  is  obtained  by  making  use  of  the  mean  temperature 
rise  profile  T(y),  the  molar  heat  capacity  of  the  carrier  gas,  Cp  ,  and 
the  heat  release  per  mole  of  reactant,  AQ .  The  product  thickness  of  the 
curves  of  Pigs.  3a,  4  normalized  by  the  1%  thickness  is  shown  in  Fig.  6a 
as  a  function  of  the  nitric  oxide  concentration.  For  small  values  of 
(NO),  there  is  a  rapid  increase  in  the  amount  of  product  with  the  kinetic 
rate  in  what  could  be  called  the  ' slow-chemistry 1  regime.  At  the  larger 
values  of  [NO],  the  product  amount  is  nearly  independent  of  [NO]  implying 
that  the  reaction  is  mixing  limited,  i.e.  that  the  chemical  rate  is  fast 
relative  to  the  mixing  rate.  Also  shown  in  Fig.  6  are  the  results  for  the 


U!  =  44  m/s  case.  The  same  data  are  shown  in  Fig.  6b  as  a  function  of  the 
Damkohler  number,  Da  (defined  as  the  ratio  of  the  mixing  time  to  the  chem¬ 
ical  time).  We  have  chosen  the  time  of  flight  at  the  mean  convection 

speed  ( U ^ -eCJ2 ) / 2  from  the  mixing  transition.  xt,  to  the  measuring  station, 
x,  to  be  an  appropriate  measure  of  the  mixing  time  (see  Broadwell  4Mungal 
[7]).  The  CHEMKIN  results  are  used  for  the  chemical  time.  This  choice 
gives  the  Damkohler  number  Da^  shown  in  Fig.  6b.  Based  on  a  transition 
Reynolds  number  of  2*10^  we  estimate  x^  =  13.6  cm  for  the  Uj  =  22  m/s  case 

and  x^  =  6.8  cm  for  the  Uj  =  44  m/s  case,  so  that  the  mixing  times  are 

20.8  msec  and  12.6  msec  respectively.  Based  on  this  criterion,  the  mixing 
limited  (fast  chemistry)  regime  occurs  past  Dax  *  40.  Finite  rate  effects 
occur  below  this  value  with  a  symmetric  profile  occuring  for  Da^  *  4 

Alternately,  one  may  choose  the  local  large  scale  turnover  time  to  be 
the  mixing  time  (ij/^u)  to  define  a  different  Damkohler  number,  Da^.  By 

this  definition,  Da^  varies  from  0.8  to  17  for  the  range  of  [NO]  investi¬ 

gated.  This  suggests  that  past  Da^  -  jq,  the  reaction  is  mixing  limited 

while  for  Da^  -  1(  the  reaction  profile  is  symmetric.  Wallace  [15]  for  a 
nitric  oxide  -  ozone  mixing  layer  was  able  to  achieve  1  <  Da^  <  16  at  high 
Reynolds  number,  while  Masutani  [16]  achieved  2  <  Da^  <  28  at  pre- 
transitional  Reynolds  numbers  using  the  same  chemical  system  as  Wallace 
A  further  discussion  of  Fig.  6b  will  be  made  in  Broadwell  4  Mungal  [7] 
when  a  detailed  comparison  of  the  results  are  made  with  the  model. 

Figures.  7a,  7b  show  sample  temperature  vs  time  traces  obtained  from 
the  eight  cold  wires  at  Uj  =  22  B/s.  for  [N0]/[N0]*  =  1  and  1/32,  i.e. 
for  fast  and  slow  chemistry.  The  time  traces  show  the  instantaneous  tem¬ 
perature  rise  recorded  by  each  probe  as  a  function  of  time,  normalized  by 
the  highest  temperature  rise  of  any  probe  during  this  time  interval 
(labeled  Tmax).  The  vertical  distance  between  a  consecutive  pair  of  hor¬ 
izontal  axes  therefore  represents  Tmax.  The  horizontal  axis  corresponds 
to  51.2  as  of  real  time.  Flow  can  be  viewed  as  being  from  right  to  left 
with  the  high-speed  fluid  at  the  top.  The  time  axis  is  greatly  compressed 
in  the  sense  that  such  a  plot  would  be  to  geometric  scale  only  if  the 


horizontal  distance  were  about  nine  times  the  distance  between  the  high¬ 
speed  and  low-speed  probes. 


We  have  discussed  such  traces  in  [1,2,4]  and  noted  the  presence  of 
large  regions  of  high  temperature,  associated  with  the  large  scale  struc¬ 
ture,  separated  by  tongues  of  unreacted,  cool  gas.  Here  however,  we  wish 
to  discuss  the  ramp-like  temperature  traces  shaded  in  Fig.  7a,  the  data 
for  the  high  kinetic  rate,  and,  to  point  out  that  they  are  not  present 
when  the  rate  is  slow,  Fig.  7b. 

If  one  pictures  a  large  structure  undergoing  reaction,  then  for  the 
present  choice  of  reactant  concentrations,  the  large  structure  is  hydrogen 
rich.  The  effect  of  free-strean  addition  of  fluid  to  the  structure  is  to 
create  a  situation  whereby  the  composition  of  the  upper  and  downstream 
edges  of  the  structure  is  biased  toward  the  high-speed  fluid,  while  the 
composition  of  the  upstream  and  lower  edges  are  biased  toward  the  low- 
speed,  lean  reactant  fluid.  This  leads  to  lateral  (y)  and  axial  (x)  gra¬ 
dients  in  composition,  and  hence  in  temperature.  In  the  fast  chemistry 
case  the  consequences  are  the  ramp-like  structures  and  time  traces  which 
on  the  lean  reactant  edge  reach  generally  higher  instantaneous  tempera¬ 
tures  than  those  of  the  rich  reactant  edge.  The  skewed  mean  temperature 
profiles  as  shown  in  Fig.  3b  result  from  long  time  averages  of  such 
traces.  These  features  are  to  be  contrasted  with  the  time  traces  of  Fig. 
7b,  for  a  slow  reaction  rate.  Here  the  ramp-like  features  are  much  less 
apparent  (become  more  top-hat)  and  the  time  traces  of  all  probes  generally 
reach  the  sane  maximum  instantaneous  temperatures  from  one  edge  of  the 
layer  to  the  other.  The  long  time  average  of  such  traces  leads  to  the 
symmetric  mean  temperature  profiles  of  Fig.  3a.  We  note  that  ramps  do  not 
occur  for  every  structure  in  the  fast-chemistry  limit  but  tend  to  be 
observed  more  often  than  not.  Similarly,  ramps  are  sometimes  present  when 
the  reaction  rate  is  low,  but  their  occurance  is  quite  infrequent. 
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We  also  note  that  since  the  fluid  mechanics  is  unaffected  by  these 
low  values  of  heat  release,  the  instantaneous  mixture  fraction  distribu¬ 
tion  of  high-speed  to  low-speed  fluid  in  the  mixing  layer  is  similar  for 
the  fast  and  slow  chemistry  cases.  In  the  fast  chemistry  limit,  the 
instantaneous  mixture  fraction  uniquely  determines  the  temperature  at 
every  point  in  the  layer  (see  Bilger  [17]  and  [2]).  For  lower  reaction 
rates,  this  is  no  longer  true  -  the  temperature  at  any  given  point  being 
determined  by  its  history.  We  also  note  that  the  two  main  features  of  the 
slow  chemistry  results  discussed  above  -  namely,  a  symmetric  mean  tempera¬ 
ture  profile  and  the  lack  of  ramp-like  structures  are  in  agreement  with 
the  measurements  of  Koochesfahani  [3]  in  a  chemically  reacting  liquid  mix¬ 
ing  layer,  in  the  limit  of  fast  chemistry,  at  comparable  Reynolds  numbers. 
The  connection  between  these  seemingly  unrelated  observations  (slow  chem¬ 
istry  in  gases  and  fast  chemistry  in  liquids)  and  a  further  discussion  of 
these  results  appear  in  Broadwell  &  Mungal  [7]. 


6.  Conclusions 

The  experiments  discussed  above  describe  the  behavior  of  the  product 
field  in  a  two-dimensional,  turbulent  mixing  layer  as  the  overall  chemical 
rate  is  varied  while  the  fluid  mechanics  and  the  heat  of  reaction  remain 
unchanged.  As  we  proceed  from  the  mixing  limited  regime  to  the  chemical 
rate  limited  regime,  three  main  effects  are  observed:  (1)  the  amount  of 
product  declines,  (2)  the  mean-temperature  profile,  which  is  initially 
somewhat  unsymmetrical  because  of  the  large  hydrogen-fluorine  reactant 
concentration  ratio,  becomes  symmetrical,  and  (3)  the  ramp-like  instan¬ 
taneous  temperature  traces  within  the  large  structures  become  more  top- 
hat.  At  the  mixing  limited  (fast  chemistry)  regime,  the  Damkohler  number 
based  on  time  of  flight  from  the  mixing  transition  exceeds  about  40. 
Below  this  value  we  enter  the  rate  limited  regime  with  a  symmetric  profile 
attained  at  a  Damkohler  number  of  about  4  (based  on  the  local  large  scale 
turnover  time,  the  Damkohler  numbers  become  10  and  1  respectively). 
Detailed  comparison  of  these  results  with  a  theoretical  model  is  presented 


in  Broadwell  &  Mungal  [7]. 
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Table 

1  -  Summary  of 

Results . 

Uj  =  22 

m/s . 

[NO] 

[NO]* 

«1 

ymax 

A 

k 

(x-x0, 

Tf  lm 

A 

*1 

k* 

2 

0.170 

0.540 

0.0457 

0.251 

1.53 

3/2 

0.163 

0.531 

0.0419 

0.231 

1.29 

1 

0.165 

0.537 

0 . 0436 

0.240 

1.00 

1/2 

0.171 

0.511 

0.0429 

0.236 

1/1.49 

1/4 

0.165 

0.488 

0 . 0408 

0.225 

1/2.37 

1/8 

0.165 

0.443 

0.0365 

0.201 

1/3.98 

1/16 

0.165 

0.396 

0.0330 

0.182 

1/6.80 

1/23 

0.169 

0.253 

0.0220 

0.121 

1/9.11 

1/32 

0.168 

0.197 

0.0169 

0.093 

1/12.2 

0 

- 

0 

0 

0 

0 

Table 

2  -  Summary  of 

Results . 

Uj  =*  44 

m/s. 

LNQ1. 

[NO]* 

«1 

ymax 

12 

k 

(x-xo) 

Tf  lm 

A 

*1 

k* 

2 

0.165 

0.491 

0 . 0405 

0.223 

1.53 

1 

0.163 

0.475 

0.0391 

0.215 

1.00 

1/4 

0.169 

0.433 

0 . 0364 

0.200 

1/2.37 

1/16 

0.168 

0.301 

0.0264 

0.145 

1/6.80 

0 

- 

0 

0 

0 

0 

Tflm  a  adiabatic  flame-temperature  rise  *  158. 7K;  Tmax  a  maximum  value  of 
T:  s  width  of  layer  where  T  =  1,1  of  Tmax;  A  »  .py jg/Td? ;  y  =  '  A 

-  nominal  condition  at  22  m/s;  (x-xQ)  =  45.7  cm  for  all  runs;  k  =  1/t  = 
overall  chemical  rate;  t  *  chemical  time  determined  from  CHEMKIN. 
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Table 

3  -  Curvefit  constants,  « 

22  m/s.  y„ 

=  -.35 

i 

■  LNQi, 

[NO]* 

co 

C1  c2 

c3 

c4 

W»™*.y-V  ,■  .;v 

2 

4.0596 

-0.9089  -0.5342 

-0.4929 

-0.6328 

3/2 

3.9804 

-0.1009  -0.4528 

-0.4400 

-0.7601 

1 

4.0568 

-0.9833  -0.5925 

-0.4336 

-0.6720 

1/2 

4.0402 

-0.9582  -0.6372 

-0.3622 

-0.5524 

1/4 

4.0497 

-0.8729  -0.6610 

-0.4596 

-0.6561 

P5 

1/8 

3.9407 

-0.9110  -0.6674 

-0.4094 

-0.6425 

1 

1/16 

3.9311 

-0.8084  -0.8073 

-0.3656 

-0.5957 

1/23 

3.5560 

-0.6492  -0.8642 

-0.4248 

-0.5362 

1/32 

3.3329 

-0.5917  -0.8614 

-0.3242 

-0.5277 

-wV 

0 

0 

0  0 

0 

0 

[£%* 

S l\< 

/y 

wtfl  '-f\  v* 

&*§■ 

Table 

4  -  Curvefit  Constants,  Uj  . 

44  n/s,  yQ 

*  -  .32 

"k\>‘ 

LNQI 

[NO]* 

Co 

C1  C2 

c3 

C4 

_ _ 

2 

4.1344 

-0.8340  -0.7880 

-0.3105 

-0.5942 

1 

4.1159 

-0.8224  -0.8642 

-0.3837 

-0.6047 

si: 

1/4 

4.0793 

-0.7320  -0.9205 

-0.2822 

-0.4830 

1/16 

3.7878 

-0.4951  -0.8463 

-0.4140 

-0.5587 

0 

0 

0  0 

0 

0 

■*, 

v* 

£y 

Mean  Temperature 

Rise,  T  (° 

K)  recorded  in  physical  coordinates,  y 

Si 

where 

& 

T(y)  -  exp 

(co  +  ciy  +  c2 y2  * 

c3y3  -  c4y 

4) 

-3  <  y  <  2.  ii  * 

(V-y0)/(x- 

xn)  and  (x-x_)  =  18.0  in.  for 

all  runs 

pcs 

il 

V7? 

\>  v 

. -■»V«  v -\V '  • 

-  •_»  “V  y  V 
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( inches ) 
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Appendix  1 


The  ideas  in  this  section  were  developed  by  J.E.  Broadwell  to  provide 
a  simple  understanding  of  how  the  NO  concentration  controls  the  overall 
reaction  rate. 


It  was  shown  above  that  reactions  ( 1 )  —  ( 3 )  can  be  considered  as  being 
representative  of  the  ^-Fg-NO  chemistry,  since  the  addit  ion  of  reactions 
( 4 ) - ( 7 )  showed  essentially  no  change  in  the  behaviour  of  the  system  as 
seen  from  Pig.  2.  It  is  possible  to  estimate  the  overall  behaviour  of  the 
system  by  the  following  arguments.  From  (2),  (3) 


dT  ‘  ‘k2»F2  '  k3N0F2 


fr  ■  -k3N°F2 


This  represents  two  equations  for  the  three  unknowns  F 2>  >J0.  H.  By  using 
either  the  rate  equations  or  conservation  of  atoms  it  is  possible  to  show 


H  +  F  +  NO  =  const  =  NO. 

l 


since  the  initial  concentrations  (subscript  i)  of  H  and  F  are  zero.  Exam¬ 
ination  of  ( 1 ) - ( 3 )  together  with  the  fact  that  would  suggest 
that  any  F  atoms  produced  (slowly)  in  (2),  (3)  would  be  (quickly)  depleted 
in  (1),  hence  a  reasonable  approximation  would  be  to  take  F  '  0.  This 
approximation  now  simplifies  (11)  to  H  >  NOj-no  which  can  be  substituted 
into  (9) ,  (10)  to  yield 


dt  *  ^2NOi^2  +  ^2-1<3^OF2  “  -1<2NOiF2  +  ^2N0F2 


V  _ 


A 


16 


=  -k3NOF2  (13) 

the  system  now  being  reduced  to  two  equations  in  the  two  unknowns  F  an(j 
NO. 


In  order  to  proceed  further,  the  CHEMKIN  simulations  were  used  as  a 
guide  as  they  suggested  two  types  of  behaviour.  In  the  first  case,  for 
saall  NO  concentrations,  CHEMKIN  suggested  that  the  NO  would  deplete  while 
the  F2  concentration  had  regained  relatively  fixed.  Thus,  if  we  assume 

-k3p2i^ 

^2  *  ^2i  we  s°lve  (13)  which  yields  NO/NOj  »  e  and  substituting 

into  (12)  gives 


dF2 

It 


-k2NOiF2(i- 


'k3F2it, 


or 


dP2  „  -k3P2i  , 
dT  "  ~F2(1"e  k2NO.r) 

k3 

where  t  »  kgNO^t.  The  second  term  can  be  neglected  for  N0t  <<  j^-5?2i 
yield  the  approximate  solution 


F2  -k2N0it 


k3 

Thus  for  NOj  <<  ^-F^,  the  overall  rate  will  be  proportional  to  NOi . 


y.y. 


The  second  case  suggested  by  CHEMKIN  concerns  large  NO  concentrations 
for  which  the  amount  of  NO  changed  slightly  but  did  not  deplete  during  the 
times  of  interest.  We  therefore  take  NO  *  NO^  in  (13).  To  solve,  we 
assume  in  (13)  that  Fg  *  F2^.  This  yields  the  solution 


further 
>i 


NO/NO,  ,  l-k0F0<t  (This  expression  would  imply  negative  NO  concentrations 


3  21 


y.t. 
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for  long  time  and  must  of  course  be  valid  for  times  t  <<  1/k  f  Alter- 

'■J  d  X 

nately  we  might  interpret  this  as  the  leading  term  in  an  expansion  for  the 
NO  behaviour).  With  this  expression.  (12)  becomes 


with  solution 


dF2 

dT  =  -k2k3NOiF2iF2t 


■2  -k2k3NOi*?2it 
—  =  e 


Thus  for  large  NO  concentrations,  the  overall  rate  is  now  proportional  to 
NO*.  The  overall  reaction  rate  derived  from  the  CHEMKIN  simulations  is 
shown  in  Fig.  8  as  a  function  of  NO  concentration.  The  initial  linear 
trend  followed  by  the  square  root  behaviour  seems  consistent  with  the 
preceeding  analysis. 

The  arguments  above  suggest  that  the  overall  rate  varies  lineariy 
with  NO,  for  small  NO,  and  then  as  the  square  root  of  NO  for  larger  values 
of  NO.  This  of  course  is  a  complication  and  raises  the  question  of 
whether  a  single  overall  rate  can  be  found  for  each  run  condition  that 
accurately  describes  the  behaviour  of  the  system.  Using  the  linear  extra¬ 
polation  technique  shown  Fig.  2  for  each  case  yields  a  characteristic 
time,  t  (see  last  column  of  Table  1).  The  CHEMKIN  results  for  all  cases 
using  reactions  ( 1 )  —  ( *7 )  are  shown  in  Fig.  9,  where  the  abscissa  has  been 
scaled  by  the  appropriate  characteristic  time.  The  similarity  of  shape  for 
all  curves  suggests  that  choosing  a  single  overall  rate  seems  reasonable 
for  describing  the  reaction  set. 
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Appendix  2 


For  the  record  we  include  the  following  curvefit  constants  for  the 
runs  reported  in  References  2,4.  The  sane  notation  as  Tables  3,4  is  used, 
and  (x-xq)  =  18.0  for  all  runs. 


Table  5  -  Curvefit  Constants  of  Reference  2,  Table  1 


* 

T 

fin 

co 

C1 

c2 

c3 

C4 

yQ 

1/8 

165 

4.2985 

-0.7807 

-0.6762 

-0.0573 

-0.5953 

-0.175 

1/4 

149 

4.2841 

-0.6418 

-0.8488 

-0.0483 

-0.5209 

-0.175 

1/2 

124 

4.2454 

-0.5327 

-0.9712 

-0.0772 

-0.4759 

-0.175 

1 

93 

4.0740 

-0.1913 

-1.0964 

-0.2538 

-0.3959 

-0.175 

1 

93 

3 . 9809 

-0.7483 

-1.2209 

-0.6152 

-0.4835 

-0.425 

2 

124 

4.3743 

-0.5202 

-1.5259 

-0.9914 

-0.5786 

-0.425 

4 

149 

4.5414 

-0.3831 

-1.6016 

-1.1379 

-0.5972 

-0.425 

a 

165 

4 . 6320 

-0.2451 

-1.7962 

-1.4960 

-0.7309 

-0.425 

Table 

6  -  Curvefit  Constants  of  Reference  2. 

Table  2 

4 

T 

fin 

co 

C1 

c2 

c3 

C4 

y0 

1/2 

137.5 

3.9154 

-1.0811 

-0.9082 

0.8426 

-0.6098 

+0 . 300 

1 

109.0 

3.8242 

-0.9398 

-1.0662 

0.5143 

-0.3058 

+0.300 

2 


154.0  4.3746  -0.7190 


-0.8895  0.3105 


-0.4022 


+0.300 
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Table  7  -  Curvefit  Contants  of  Reference  4 


Run 

c 

0 

C1 

c2 

c3 

c4 

y0 

RE1 

4.1019 

-0.7526 

-0.5933 

-0.5687 

-0.6048 

-0.325 

Re2 

4 . 1603 

-0.9880 

-0.9188 

-0.5182 

-0.6362 

-0.325 

Re3 

4.1105 

-0.7930 

-0.9572 

-0.6126 

-0.6033 

-0.325 

Re2T 

4.0359 

-1.7096 

-1.1633 

-1 . 0232 

-3.4653 

-0.210 

R63T 

4.0223 

-1.7280 

-2.0311 

-0.7710 

-1.9014 

-0.210 

R64 

3.9053 

-1.7094 

-1.5894 

-0.8499 

-1.9001 

-0.210 

y, y. y.y.y.y v.-. 


Sa& 


-*>  s 
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9.  Figure  Captions 


Fig  Turbulent  Shear  Layer  Geometry. 

Fig  2  CHEMKIN  Simulation  of  Constant  ?~essure  Reactor  using  Reactions 
{ 1 )  -  ( 3 )  and  (i)-(T).  [NO] /[NO’*  =  1. 

rig.  3a  Normalized  Mean  Temperature  Profiles  for  decreasing  .’NO., 

22  m/s.  [NO] /[NO]  •  □  -  1 .  O'  1/2.  A  *  1/4  ■+■  -  i  3 ,  X  - 

1/16,  V  -  1/23.  O  -  1/32,  ♦  -  C. 

Fig.  3b  Scaled  Mean  Temperature  Profiles  for  decreasing  [NO],  U,  -  22 

m/s.  Symbols  same  as  Fig.  3a. 

* 

Fig.  3c  Scaled  Mean  Temperature  Profiles  for  [NO]/[NO]  =  1  &  1,32,  ’J .  - 

22  m/s.  Symbols  same  as  Fig.  3a. 

Fig.  4  Normalized  Mean  Temperature  Profiles  for  increasing  [NO],  U  = 
22  m/s.  [NO]/ [NO j  :  □  -  1.  O  ~  3/2,  A  -  2.  1 

Fig.  5a  Normalized  Meajp  Temperature  Profiles  for  varying  [NO],  U  =  u 

m/s.  [NO]/ [NO]  :  □  -  1,  O  *  1/4.  A  "  1'16.  +  -2.  1 

Fig.  5b  Scaled  Mean  Temperature  Profiles  for  [NO]/[NOj*  -  1  &  1/16,  F. 

44  m/s.  Symbols  same  as  Fig.  5a. 

Fig.  6a  Normalized  Prcduct  Thickness  vs  Nitric  Oxide  Concentration. 

Fig.  6b  Normalized  Product  Thickness  vs  Time  of  Flight  Damkonier  Number. 

Fig.  7a  Temperature  vs  Time  Trace,  [NO] /[NO]  -  1.  U,  -  22  ms. 

Fig  7b  Temperature  vs  Time  Trace,  [NO].'[NC]  =  1/32,  C,  =  22  m/s. 

Fig  8  k/k  vs  [NO]/[NO]  from  CHEMKIN  simulations  for  Reactions  ; 

(7)  . 

Fig  9 


CHEMKIN  Simulations  of  Constant.  Pressure  Reactor  using  Reactions 
( 1 ) — ( 7 )  for  [NO]  shown  in  Table  1. 
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Abstract 


A  simple  model  describing  chemical  reaction  in  a  turbulent  shear 
layer  is  put  in  quantitative  form  and  the  model  predictions  compared 
with  experiment.  The  reactants  are  not  pre-mixed  and  the  flow  is 
two-dimensional.  The  dependence  of  the  amount  of  product  in  the  layer 
on  the  Schmidt,  Reynolds,  and  Damktthler  numbers  and  on  the 
stoichiometric  ratio  is  exnibited  explicitly  in  the  model.  The  first 
two  parameters  appear  as  the  product  (Sc  Re)  and  influence  that  part  of 
the  reaction  taking  place  in  strained  laminar  flames.  In  the  limit 
Sc  Re  +  «  (for  Sc  >  1  and  Re  »  1)  the  molecular  mixing  becomes 
independent  of  these  parameters  as  it  is  in  classical  turbulence  theory 
but  the  composition  and  spatial  distribution  of  the  molecularly  mixed 
fluid  does  not  resemble  those  given  by  such  theories.  The  model 
predictions  are  in  reasonable  agreement  with  the  experimental  results 
for  a  DamkOhler  number  range  from  zero  to  a  mixing  limited  value.  The 
comparison  with  a  set  of  data  for  gases  in  which  Re  was  varied  is 
somewhat  less  satisfactory. 
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I.  Introduction 


This  paper  presents  the  results  of  an  analytical  study  of  chemical 
reactions  in  a  two-dimensional  turbulent  shear  layer.  In  particular, 
the  experimental  results  for  the  variable  rate  hydrogen-fluorine 
reaction  reported  by  Mungal  and  Frieler  in  Ref.  1  are  interpreted  in 
terms  of  the  model  proposed  in  Ref.  2.  The  reaction  takes  place  between 
H2  and  nitric  oxide,  NO,  carried  in  a  high  speed  nitrogen  stream  and  F2 
in  a  low  speed  nitrogen  stream.  In  most  of  the  experiments  the  high  and 
low  speed  velocities  were  22  m/3ec  and  8.8  m/sec  but  several  runs  were 
made  at  44  m/sec  and  17.6  m/sec.  The  H2  concentration  is  eight  percent 
molar,  the  F2  concentration  is  one  percent  molar,  and  that  of  the  NO 
(needed  to  initiate  the  reaction)  varies  between  0  and  0.03  percent 
molar.  A  description  of  the  apparatus  and  of  the  experiment  is  given  in 
Ref.  1. 

The  important  reactions  are: 


ki 

F  +  H2  •»  HF  ♦  H 
k2 

F2  +  H  **  HF  +  F 

ko 

F2  +  NO  -o'3  NOF  +  F 

describet  by  the  equations: 


(1) 

(2) 

(3) 
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-  K,  [F]  [H2] 


d[H2] 
dt 

dCF2] 

— — —  ■  ~ k2  [F2]  [H]  ~  [NO]  [F2] 


d[NO] 

dt 


-k3  [F2]  [NO] 


It  was  shown  in  Ref.  1  that  this  kinetic  system  behaved  as  if  it 
were  a  single  reaction 


k 

F2  -  2 


HF 


(4) 


governed  by  the  rate  equation 


d[HF  ] 
dt 


k  [H2][F2] 


(5) 


where  k  is  a  function  of  the  initial  nitric  oxide  concentration,  [N0j_]. 
In  the  following  we  make  use  of  this  approximation,  i.e.,  we  assume  a 
single  reaction  governed  by  a  rate  constant  k. 

Because  the  reactant  concentrations  are  low,  the  density  changes 
are  neglected  and  the  fluid  motion  is  assumed  to  be  unaffected  by  the 
chemical  reaction. 
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II.  Physical  Ideas 


The  ideas  on  which  the  model  is  based  come  from  a  series  of 
experiments  beginning  with  those  of  Brown  &  Roshko’.  Most,  in  fact, 
derive  from  this  now  widely  known  seminal  study.  The  subsequent  work  of 
Konrad*,  Breidenthal*  and  Brown  &  Dimotakis*  and  the  analyses  of  Corcos 
and  Sherman’  provided  further  evidence  for  the  basic  concepts.  This 
evidence  is  discussed  in  detail  in  Ref.  2  and  a  model  suggested  by  the 
observations  is  proposed  in  qualitative  terms.  Here  the  physical  ideas 
will  be  reviewed  briefly,  the  model  will  be  put  in  quantitative  form, 
and  results  from  it  compared  with  the  experimental  findings.  Finally, 
shortcomings  of  the  model  are  discussed  and  possible  improvements 
suggested. 

The  experiments  cited  above  show  that  the  sequence  of  events  by 
which  parcels  of  fluid  from  the  two  streams  reach  the  molecularly  mixed 
state  begin  when  the  fluid  elements  are  entrained  into  the  shear  layer 
by  large  scale  coherent  motions.  The  parcels  are  distributed  across  the 
layer  and  are  entangled  (using  Roshko's  term)  by  vortical  motions  in 
large  scale  structures.  This  motion,  molecular  mixing,  and  chemical 
reaction  are  governed  by  the  equations,  in  non-dimensional  form: 


+  V*  •  7C4 


(8) 


P 

£, 


a 


48 


V2ai  -  co* ^ 


in  which  V*  -  V/Au,  p*  -  p/p„(Au)2,  x^  -  x^s,  t*  -  t(A)/5),  Sc  -  v/D, 
Re  -  Au  5/v,  ai  is  the  mole  fraction  and  <^*  the  dimensionless  chemical 
reaction  term.  Here  Au  is  the  velocity  difference  (u1  -  113)  and  5  the 
local  layer  thickness.  We  take  the  Schmidt  number,  Sc,  to  be  of  order 
unity  or  larger  and  the  Reynolds  number,  Re,  to  be  large.  It  is 
conventional,  and  presumably  justifiable,  at  large  Re  to  ignore  the 
effects  of  the  viscous  terms  in  Eq.  (7)  on  the  motion.  For  the 

molecular  mixing,  however,  the  diffusion  term  is  essential.  In  view  of 
the  nature  of  Eq.  (8),  especially  in  the  limit  Re  »,  it  seems 
plausible  that  thin  diffusive  layers,  boundary  layers  in  the 

mathematical  sense,  would  form  between  fluid  from  the  two  streams  and 
that  their  thickness,  t,  would  be  of  the  order  6/(ScRe)1/2.  In  the 

model,  these  layers  are  represented  by  the  solutions  for  the  laminar 

diffusion  zones  formed  in  stagnation  point  flows,  conventionally  called 
strained  flames. 


It  is  argued  next  that  as  a  structure  moves  downstream,  the  laminar 
mixing  layers  are  lengthened  and  further  intertwined  as  ever  smaller 
scales  are  generated.  This  (Lagrangian)  motion  is  taken  to  be  a 
Kolmogorov-like  cascade  to  Xq,  the  Kolmogorov  scale.  When  the 
fluctuations  reach  this  scale,  the  layers  merge.  (During  this  motion, 
the  structures  are  distorted  and  may  intermingle,  so  that  when  the 
layers  merge  they  do  not  necessarily  form  large,  distinguishable  regions 
in  the  cores  as  Fig.  5  in  Ref.  2  may  unfortunately  suggest.)  Since  the 
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is  of  the  order 


time  to  diffuse  across  the  scale  Xg  -  6/Re^/4. 

(6/Au)  •  Sc/Re1 /^l  at  high  Re  the  molecules  in  the  layers  interdiffuse  in 
a  time  that  is  short  compared  to  the  cascade  time,  6/Au,  to  form  what  is 
called  in  the  model,  the  homogeneous  mixture.  Because  the  layers  in  a 
given  structure  have  been  distributed  throughout  the  mixing  layer  by  the 
large  scale  motion  before  they  merge,  we  assume  that  the  homogeneous 
mixture  has  a  composition  that  is  nearly  independent  of  the  lateral 
coordinate,  y.  Additional  experimental  evidence  supporting  this 
approximation  is  contained  in  the  results  of  Mungal  &  Dimotakis*,  and 
especially,  Koochesfahani *  and  Koochesfahani  &  Dimotakis10.  In  the 
model,  as  a  further  idealization,  the  homogeneous  mixture  composition  is 
taken  to  be  a  single  value:  it  is  assumed  to  be  composed  of  fluid  from 
the  two  streams  in  the  ratio,  E,  the  ratio  of  high  speed  to  low  speed 
fluid  entrained  into  the  mixing  layer. 

The  flame  sheet  area  per  unit  volume,  S,  computed  in  the  cascade  by 
taking  S  -  1  / X ,  where  X  in  the  fluctuation  scale,  for  6  >  X  >  Xq, 
behaves  as  intuition  suggests,  i.e.,  S  increases  slowly  at  first  and 
then  extremely  rapidly  as  X  +  Xq.  This  behavior  is  the  motivation  for 
dividing  the  molecularly  mixed  fluid  into  two  distinct  states,  flame 
sheets  and  the  homogeneous  mixture.  In  the  model,  then,  S  is  the  flame 
sheet  area  per  unit  volume  at  the  large  scale  and  is  taken, 
consequently,  to  scale  with  1/6  and  to  be  independent  of  Re.  (The 
thickness,,  however,  as  already  noted,  depends  on  Re  as  well  as  Sc.) 
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This  description  of  the  flow,  involving  as  it  does  the  motion  of 
large  scale  fluid  elements,  pictures  the  flow  to  be  unsteady  at  spatial 
and  temporal  scales  6  and  6/Au.  Evidence  for  such  unsteadiness  is  clear 
in  the  temperature-time  traces  in  Ref.  1.  In  the  further  development  of 
the  model,  however,  we  will  be  concerned  only  with  mean  values  and  will 
put  the  ideas,  which  were  most  naturally  described  in  Lagrangian  terms, 
into  Eulerian  form. 

The  ideas  presented  above  can  be  summarized  and  clarified  with  the 
help  of  the  diagram  in  Fig.  1  which  is  intended  to  represent  the  three 
fluid  states  in  the  shear  layer  and  to  show  the  sequence  by  which  fluid 
elements  pass  from  one  to  the  other.  Also  indicated  explicitly  is  the 
presence  of  unmixed  fluid  within  the  layer.  The  v's  are  volume  fluxes 
per  unit  length  into  the  various  states,  with  v1  and  v2  the  fluxes  from 
the  high  and  low  speed  streams  into  the  layer.  (The  entrainment  ratio, 
E  (>  D,  is  v1/v2.) 

When  the  molecular  diffusion  coefficients  of  the  contituents  are 
equal,  the  flame  sheets  consist  of  equal  volumes  from  the  two  streams. 
The  excess  high  speed  fluid,  therefore,  becomes  part  of  the  homogeneous 
mixture,  directly,  at  the  end  of  the  cascade. 

The  homogeneous  mixture  volume  flux  increases  with  axial  distance  x 
by  the  addition  of  a  volume  flux  per  unit  length,  vh,  the  sum  of  the 
flame  sheet  flux,  vf  and  of  high  speed  fluid,  vd.  It  is  easy  to  show 


The  justifications  for  the  assumption  that  vh  is  a  constant 


independent  of  Re  (to  be  determined,  implicitly,  by  experiment  as  shown 


later)  are  discussed  in  Ref.  2;  these  are,  essentially,  the  same  as 


those  given  above  for  dividing  the  mixed  fluid  into  the  two  states.  The 


result  is  that  the  amount  of  homogeneously  mixed  fluid  rises  linearly 


with  x  as  does  the  entrained  fluid. 


This  description  of  the  flow  is  shown  schematically  in  Fig.  2  in 


which  the  cross-hatched  regions  represent  the  homogeneous  mixture  and 


the  thin  lines  the  flame  sheets.  In  the  present  formulation  it  is 


assumed  that  the  flame  sheets  are  formed  only  between  pure  fluid  from 


the  two  streams  and,  in  addition,  that  their  structure  and  composition 


is  in  quasi-equilibrium  with  the  local  large  scale  strain.  (In  Ref.  2 


it  is  wrongly  stated  that  these  properties  are  in  equilibrium  with  the 


strain  at  the  small  scales.) 
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III.  Model  Formulation 


In  the  experiments  to  be  discussed,  the  amount  of  product  in  the 
shear  layer  is  characterized  by  the  "product  thickness",  6p,  defined  at 
a  particular  axial  measuring  station  by 


tm  Cp(y) 

J  "E“dy 


where  Cp  is  the  molar  product  concentration  and  C*  i3  the  reactant 
concentration  in  one  of  the  streams,  here  taken  to  be  the  low  speed 
fluorine  bearing  stream.  Following  the  discussion  above,  the 
contributions  of  the  flame  sheets  and  the  homogeneous  mixture  are 
considered  separately  as  follows: 


r“  (cp)n  (Cp)f 

11  r o(x,y)  dy  +  — —  l  •  S  dy  (9) 

_  _ 


where  (Cp)h  and  (Cp)f  denote  the  homogeneous  mixture  and  average  flame 
sheet  product  concentrations,  in  moles  per  unit  volume  of  homogeneously 
mixed  fluid  and  of  flame  sheet  fluid ,  respectively.  The  fractional 
volume  of  homogeneously  mixed  fluid  is  denoted  by  a(x,y)  and  that  of  the 
flame  sheet  by  l  •  S  where  l  is  the  sheet  thickness. 


noted,  (Cp)h  is  independent  of  y  at  fixed  x  and 
a(x,y)  -  a(n).  where  n  •  yA5  and  <5  is  the  shear  layer  thickness,  here 
taken  to  be  distance  between  the  one  percent  of  maximum  of  the  product 
profiles.  The  homogeneous  mixture  term  can,  consequently,  be  written 
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Carrier,  Fendell,  and  Marble11  show  that  l  -  a(D where  a  is  a 
constant  of  order  one,  D  is  the  molecular  diffusion  coefficient,  and  e 
the  strain  rate.  As  stated  above,  we  take  e  to  be  the  large  scale 
strain,  Au/5  and  assume  that  l  and  (Cp)^  are  also  independent  of  y  at 
fixed  x.  Since  S  scales  with  1/6,  we  put  S  -  B(n)/6  and  write  for  the 
flame  s'  -at  term 


(Cp)f 

C_ 


»)f  f"  (Cp)#  .nx,1/2  fm 

—  I  j  8(n)  dn  -  a  — - —  ( ~)  j  B(n) 


a  •  (Cp)f 


S(ScRe)“1/2  |  B(n)  dn 


With  these  expressions,  Eq.  (9)  becomes 


y  -  A  (Cp)h/C.  ♦  B  (Sc  Re)~1/2  (CP)f/C. 


cm  m 

where  A  -  [  a(n)  dn  and  B  -  a  B(n)  dn  are  constants  to  be 


determined  as  described  later. 


Observe  that  A  and  B  have  to  do  only  with  the  molecular  mixing; 
all  the  chemical  kinetic  effects  appear  in  (Cp)h  and  (Cp)f.  Strictly,  A 
and  B  must  be  determined  for  each  free  stream  velocity  and  density 
ratio.  It  is  likely,  however,  from  the  arguments  leading  to  the  model, 
that  they  are  weak  functions  of  these  variables. 
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In  the  next  two  sections,  methods  for  determining  (Cp)h  and  (Cp)f 
are  derived,  with  the  flame  sheets  considered  first. 
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IV.  Strained  Flames 


A  two-dimensional  strained  flame  is  formed  when  two  streams 
carrying  reactants  meet  at  a  stagnation  point  as  sketched  in  Fig.  3. 
Under  our  conditions,  the  overall  conservation  equation  is 


3u  +  3v 
3x  3y 


0 


If  v  is  a  linear  function  of  y. 


v--cy  .  fs.-e  <"> 

When  e  is  constant  and  the  reaction  is  fast  relative  to  diffusion,  the 
reaction  takes  place  on  a  plane  parallel  to  x  and  the  exact  solution  of 
Ref.  11  is  available.  For  finite  reaction  rates  numerical  computation 
is  required1*.  While  such  computations  could  be  used  here,  it  is 
instructive,  instead,  to  use  the  approximate  solution  developed  next. 

A  boundary  layer  Integral  approach  applied  to  the  control  volume 
indicated  by  dashed  lines  in  Fig.  3  yields  the  desired  results.  To 
obtain  a  simple  approximate  solution,  we  specialize  from  the  beginning 
to  large  values  of  the  ratio  (CH)„/(cF).  -  1/0,  where  0  is  the 
stoichiometric  ratio  (following  the  conventions  of  Ref.  8)  and  where  the 
subscripts  H  and  F  denote  hydrogen  and  fluorine.  Further,  consider  only 
average  values  for  the  concentrations  in  the  control  volume.  Mote  first 
that  from  Eqs.  (11)  v(-  i/2)  -  e  i/2  and  u  -  e  A  x. 
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The  conservation  equations  for  H2  and  F2  in  the  control  volume  are 


-  v(  + 1/2)  (CH)„  Ax  -  k(Cp)f  (CH)f  •  £  •  Ax  -  (CH)f  e  £  •  Ax  (12) 
v( -  1/2)  (Cp).  Ax  -  k(Cp)f  (CH)f  £  •  Ax  -  (CF)f  e  £  •  Ax  03) 


From  Eq.  (12), 


(CH>f 


1 

2  [1  ♦  *  <Cp)f ] 


Putting  this  expression  into  Eq.  (13),  and  solving  for  (Cp)^  in  the 
limit  (Cp)„  -►  0  (to  approximate  (CH)»/(Cp),  »  1)  we  get 


<CP) 


F'f 


<Cp). 


[1  . 

L  2  e  J 


In  the  same  limit. 


(Cu) 


f 


<c»). 


(14) 


(15) 


The  product  concentration,  (Cp)f,  is  obtained  from  the  control 
volume  concentration  equation 

**^p)f  (Cfl)f  Ax  •  £  »  u(Cp)f  £  *  e  Ax(Cp)f  £ 

which  becomes,  using  Eq.  (14)  and  (15), 
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k(CH) 


(16) 


(Cp) 


(Cp) 


F 


p J  r  " 


2  e 


a  result  consistent  with  the  requirement  that 

(Cp)^  *  (Cp) f  -  (Cp)./2  . 

To  compare  these  results  with  the  exact  solution  of  Ref.  11,  note 
first  that  since  l  »  a(D/e)^2,  the  amount  of  product  per  unit  flame 
sheet  area,  kp  -  (Cp)f  •  i  is,  from  Eq.  (16)  with  k  ♦  •,  given  by 


•cp-  (a/2)  (CF).  (D/e) 


1/2 


(17) 


The  exact  solution  is 


«p  -  (2/ir)1/2  (♦♦I)  e"*2  (Cp).  (D/e)1/2 


(18) 


where  A  -  erf-1  ^  y;-  .  Thus  the  dependence  on  the  parameter  (D/e)  is 

\  I  ) 

the  same  in  the  two  expressions. 


For  small  k,  Eq.  (16)  reduces  to 


(Cp) 


(Cjj).  (Cp). 


P'f 


(k/e) 


In  summary,  Eqs.  (15)  and  (16)  are  approximations  for  all  k  for 


(CH)„/(CF).  -  1/#  »  1,  while  Eq.  (18)  is  the  exact  solution  for 
k(CH). 

— - —  »  1  for  all  p. 

Arguments  were  given  above  for  taking  e  proportional  to  Au/6. 
Variation  of  the  proportionality  constant  around  unity  had  so  little 
effect  on  the  overall  results  that  unity  was  chosen  for  convenience  in 
the  following  work. 
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V.  The  Homogeneous  Mixture 

In  this  section  expressions  for  the  reactant  and  product 
concentrations  in  the  homogeneous  mixture  are  developed.  First,  an 
expression  for  the  flux  of  the  homogeneously  mixed  fluid  is  written  and 
then  the  fluorine  conservation  equation  is  derived. 

The  homogeneous  mixture  flux,  Vh,  can  be  written 

m  m 

Vh  -  I  a(y)  u(y)  dy  -  6  at(n)  u(n)  dn 

—  m  —m 

where  a  has  the  same  meaning  as  before.  To  avoid  introducing  an 
additional  unimportant  constant,  we  approximate  the  above  integral  by 
writing 

OD 

Vh  -  5  u  j  Ci(n)  dn  *>  6  u  A  (19) 

where  u  -  (u-|  +  u2)/2.  Therefore,  this  flux  increases  linearly  with  x 
and  is  equal  to  vhx  as  indicated  In  Fig.  1,  i.e., 

<5  u  A  -  vh  x  (20) 

When  fluorine  is  carried  in  the  low  speed  stream,  no  fluorine  "enters" 
the  flux  directly,  it  Is  homogenized  only  after  becoming  part  of  a 
flame  sheet.  The  fluorine  and  hydrogen  conservation  equations  for  the 
homogeneous  mixture  are,  under  these  conditions, 
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—  |  (CF)h  a(y)  u(y)  dy  -  -  k(CH)h  (CF)h  a(y)  dy 


(21  ) 


2(Cf)^> 

"FT  vh 


_d_ 

dx 


a) 

(CH)h  u(y)  dy  -  -  |  *<(CH)h  (CF)h  a(y)  dy 


(22) 


*■  t  e+i  ^cH^f  +  (e+1)  vh 


Equation  (21)  deacribea  the  change  in  the  fluorine  flux  in  the  mixture 
cauaed  by  the  conaumption  of  ?2  in  the  chemical  reaction  and  the 
addition  from  the  flame  sheets.  Using  Eq.  (19)  and  (20),  we  find, 


_q_  -k(CH)h  (Cp)^  x  2(CF)f 

dx  [x^cF^hJ  “  " 


E+1 


(23) 


to  fx(CH}h] 


(24) 


-  k<CH)h  (CF)h  x 


I  TeTTT  (cH}f  +  tIttt  (ch)-1 


The  approximation  used  in  the  flame  sheet  analysis,  i.e.,  that  the 
hydrogen  depletion  by  chemical  reaction  can  be  neglected  is  introducted 
next.  With  (CH)F  -  (C^)m/2,  Eq.  (24)  yields 


) 


H  'h 


(E*1)  H 


(Cu)  , 


(25) 


The  hydrogen  from  the  free  stream  is  again  simply  diluted  in  the 
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homogeneous  mixture  by  the  low  speed  free  stream  fluid. 
With  this  result  Eq.  (23)  becomes 


[<x*(CF)h] 


E 

(E*1 ) 


2(Cp)f 

(E+1) 


(26) 


where  x*  -  k  (CH)  x/u. 

n  a* 

Since  e  -  Au/6,  in  which  Au  is  simply  related  to  u  and  5  -  0.l6x, 
Eq.  (Ill)  shows  that  (C p)f,  also,  is  a  function  of  x*  only. 
Equation  (26)  can,  therefore,  be  solved  by  standard  methods  to  yield 
( Cp ) h  aa  a  function  of  x*  and,  since 


(CpVc-  +  iCpJ./C.  -  1/(E  ♦  1)  , 


(27) 


the  variation  of  the  product  concentration  at  the  (fixed  x)  measuring 
station  with  the  reaction  rate.  Interestingly,  the  over-all  solution  is 
dependent  only  upon  the  DamkOhler-like  variable,  x*,  which  arises 
naturally  in  the  analysis. 

It  is  helpful  to  note  that  when  x/u  is  replaced  by  t  in  Eqs.  (23) 
and  (24),  they  become  equations  describing  exactly  the  fluorine  and 
hydrogen  concentrations  in  a  homogeneous  (ie.,  well-stirred)  reactor  to 
which  streams  of  flame  sheet  and  free  stream  fluid  are  being  steadily 
added. 


19 


4  June  1986 


VI.  Comparison  of  Theory  &  Experiment 


To  compare  the  theoretical  predictions  with  the  experimental 
results  of  Ref.  1,  values  for  the  constants  A  and  B  are  required.  They 
have  already  been  determined  by  Mungal  &  Dimotakis*  from  their  data  for 
fast  chemical  reactions.  For  completeness,  the  constants  will  be 
rederived  here,  in  a  slightly  different  way,  with  essentially  the  same 
results.  In  any  even*:,  no  re-adjustments  are  made  (nor  are  possible)  to 
fit  the  variable  chemical  rate  data  of  Ref.  1. 


The  most  straight-f orward  and  accurate  procedure  makes  use  of  two 
sets  of  results:  (1)  the  first  from  a  run  for  a  nominal  NO 

concentration  in  Ref.  1  and  (2)  from  the  study  of  Koochesfahani  & 

Dimotakis1®  of  a  fast  chemical  reaction  in  water.  We  need  not  be 

concerned  here  with  the  details  of  the  later  experiment,  needing  only  to 

note  that  it  was  for  Re  -  2.2  •  10**,  $  -  1/10,  and  Sc  -  600. 


The  HF  experiment,  for  Re  *  6.6  •  101*  and  <p  -  1/8  yielded 
<5 p/ 5  -  0.2U  and  the  water  experiment  gave  5p/5  -  0-13- 


Eq.  (10)  repeated  from  Section  III,  i.e., 


A  (Cp)h/C.  ♦  B(ScRe)'1/2  (Cp)f/C„ 


(10) 


is  applicable  to  both  experiments. 


20 


M  June  1986 


For  the  liquid  reaction  k  »  1  and  Eqs.  (26)  and  ( 1 6 )  can  be  used 


to  find  that  (Cp)h/C.  »  1  / ( E  +  1)  and  (CP)f/C.  -  0.5.  For  the  F2/H2 
fast  reaction  experiment,  kCC^)*  =*  (0.53  •  10"^  sec)-1  and  t  *  176  sec-1. 
We  have,  again  using  Eqs.  (26),  (27)  and  (16),  (Cp)h/c.  -  .997/(E  +  1) 
and  (Cp)f/cm  -  0.42.  With  these  values  and  E  -  1.3,  Eq.  (10)  applied  to 
the  two  experiments  yields  A  -  0.28  and  B  =  64.5. 

Substitution  of  these  values  into  Eq.  (10)  shows  that  the  product 
in  the  gas  at  this  Reynolds  number  is  almost  equally  divided  between  the 
homogeneous  mixture  and  the  flame  sheets. 

The  expression  for  finite  chemical  reaction  rate  is,  therefore, 

Y  *  0-28  (CP)h/C.  ♦  64.5(Re)“  (Cp^/C. 

where  Sc  has  been  put  equal  to  unity  and  the  quantities  (CpJ^/C*  and 
(Cp)f/C.  au-e  to  be  found  from  Eqs.  (26),  (27)  and  (16).  Figure  4 
compares  the  analytical  and  experimental  results.  It  shows  the 
normalized  product  thickness  at  the  measuring  station  in  terms  of  the 
DamkOhler  number  Da,  where  Da  »  k  s/u  and  s  is  the  distance  between 
the  measuring  station,  x  -  47.7  cm,  and  the  location  of  the  shear  layer 
transition  to  fully  developed  turbulence,  x  -  15  cm.  The  contributions 
of  the  homogeneous  zones  and  the  flame  sheets  indicated  separately,  show 
that  in  contrast  to  the  homogeneous  mixture,  the  product  concentration 
in  the  flame  sheet  approaches  its  asymptotic  value  slowly:  even  at 
Da  -  60,  the  reaction  is  not  quite  mixing-limited.  While,  as  can  be 
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seen,  the  agreement  between  the  experiment  and  the  predictions  of  the 
model  is  satisfactory,  experiments  over  a  wider  range  of  conditions  are 
needed  to  test  the  model  more  stringently. 

It  is  useful  next  to  discuss  effects  of  the  Reynolds  number  on  this 
reaction  because  when  Re  is  increased  by  increasing  the  velocities, 
there  is  an  accompanying  change  in  the  DamkOhler  number.  Fig.  5  shows 
the  theoretical  results  for  several  values  of  the  reaction  rate  together 
with  the  recent  measurements  of  Mungal  et  al.11*  The  curve  marked 
k/k*  »  1  is  to  be  compared  with  the  experimental  values.  That  for 
infinite  reaction  rate,  k/k*  -  ®,  illustrates  the  effect  of  Reynolds 
number  alone.  The  Reynolds  number  range  in  the  experiments  is  too 
limited  to  permit  definite  conclusions  from  the  comparison  at  this  time. 

The  analysis  in  this  paper  has  been  concerned  primarily  with  the 
effects  of  DamkShler  and  Reynolds  number  on  chemical  reactions  in 
turbulent  shear  layers.  The  application  of  the  model  to  reactions  in 
which  the  equivalence  ratio  was  varied  is  discussed  by  Mungal  and 
Dimotakis. * 
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VII.  Concluding  Remarks 


We  conclude  with  a  tentative  explanation  of  the  data  in  Fig.  7  of 
Ref.  1  showing  that  the  slope  of  the  ramps  in  the  temperature-time 
traces  within  the  structures  at  high  chemical  reaction  rate  decreases  as 
the  rate  is  lowered  becoming  nearly  zero  at  the  lowest  rate. 
Simultaneously,  the  asymmetry  in  the  time-mean  temperature  profile 
disappears:  Fig.  3a,  Ref.  1.  Observe  first  from  Fig.  4  that  the 
relative  contribution  of  the  flame  sheets  to  the  product  declines  as  the 
reaction  rate  is  reduced.  Next  Koochesfahani 1 1  found  that  in  liquids, 
where  (in  the  present  view)  the  flame  sheet  contribution  is  negligible, 
the  time-mean  profiles  of  the  product  concentration  are  symmetric  and 
ramps  are  not  present  in  the  concentration-time  traces.  These 
observations  are  consistent  with  the  idea  that  the  ramps  and  asymmetries 
are  associated  with  flame  sheets. 


In  the  model,  the  sheets  are  formed  between  pure  fluids  from  the 
two  streams,  consequently  their  composition  is  the  same  everywhere. 
This  over-simplif ication  must,  therefore,  be  dropped  if  the  observations 
described  above  are  to  be  explained  by  reference  to  the  flame  sheet 
composition.  The  required  generalization  is  that  the  newly  entrained 
fluids  form  mixing  layers  not  only  between  themselves  but  with  the 
homogeneous  fluid.  Mungal  and  Dimotakis*  suggest  such  a  generalization. 
After  reviewing  the  evidence  that  high  speed  fluid  enters  the  structures 
along  their  high  speed  and  downstream  boundaries  and  the  low  speed  fluid 
along  the  low  speed  and  upstream  edges,  they  propose  that  the  two 
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streams  mix,  to  some  degree,  with  the  homogeneous  fluid  as  they  enter. 
The  resulting  flame  sheet  compositions  are  such  as  to  explain  the 
observed  ramp  and  mean  profile  characteristics. 

These  ideas  may  be  stated  another  way:  at  finite  Reynolds  number, 
mean  product  profiles  are  symmetric  and  there  are  no  ramps  when  the 
molecular  mixing  is  delayed,  i.e.,  does  not  occur  in  the  flame  sheets. 
In  gases,  this  can  be  the  consequence  of  a  low  reaction  rate;  in 
liquids,  slow  molecular  diffusion  is  the  cause.  In  the  limit 
(Sc  Re)  •*  •  (with  Sc  >  1  and  Re  »  1)  the  flame  sheet  thickness  goes  to 
zero;  consequently  liquids  and  gases  behave  similarly.  That  the  model 
approaches  a  limit  in  which  the  mixing  is  independent  of  Re  and  Sc  is  in 
accord  with  classical  turbulence  theory  but  the  distribution  and 
composition  of  the  mixed  fluid  are  qualitatively  different. 

Ramp-like  features  have  been  observed  in  scalar  time  traces  in  a 
wide  variety  of  experiments  in  the  laboratory 1  ** 1  *  and  in  the 
atmosphere17.  A  more  extensive  review  of  these  earlier  studies  and  a 
comparison  with  the  present  work  is  in  preparation. 
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Figure  2.  Strained  flame  sheets  (solid  lines)  and  homogeneously 
mixed  zones  (cross-hatched)  in  shear  layer. 


29 


